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The best performing finned tube had a fin spacing of l. 
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It showed an enhancement over the smooth tube of between 4 
to 5 under vacuum conditions and around 6 at atmospheric 
conditions. Thus, the use of finned tubes may result ina 
Significant reduction in the size and weight of present-day 


steam condensers. 


tl 


i T.. 


IV. 


MON (Cincy ° Urbis wi 


ddadd + Wu's 


E OF CONTENTS 


TABL 


INTRODUCTION ------------------------ ------------ 
A. FILMWISE CONDENSATION ON FINNED TUBES ------- 
B. SURVEY OF ANALYTICAL MODELS AND 
EXPERIMENTAL WORK --------------------------- 
C. OBJECTIVES -------------------------------- == 
DESCRIPTION OF APPARATUS ------------------------ 
A. SYSTEM OVERVIEW ----------------------------- 
B. SYSTEM INSTRUMENTATION ---------------------- 
C. DATA ACQUISITION SYSTEM --------------------- 
D. SYSTEM MODIFICATIONS ----9--------9-- 
E. VACUUM INTEGRITY ---------------------------- 
F. TUBES TESTED -------------------------------- 
G. INSERT -------------------------------------- 
SYSTEM OPERATION AND PROBLEMS ------------------- 
dee ROL won CONDE Nonl LON PROBLEM =----------— 
B. VAPOR VELOCITY LIMITATIONS ------------------ 
DATA REDUCTION ----------------------------- ----- 
A. PROGRAM SIEDER ------------------------------ 
B. MODIFIED WILSON PLOT PROGRAM ---------------- 
Ca MATIN BATA REDUCTION PROGRAM (DRP4) ---------- 
RESULTS AND DISCUSSION -------------------------- 
A. INSIDE HEAT-TRANSFER COEFFICIENT ------------ 
lL. Instrumented Tube Results --------------- 
2.  Uninstrumented Tube Results ------------- 
cee Summary Of tne Inside Heat-Transfer 


Coefficient 


yi 


18 


Ze 
28 
30 
50 
Se 
34 
35 
38 
39 
42 
45 
47 
49 
uk 
53 
Shs 
56 
a 
a 
De 


64 


66 


B. PERFORMANCE OF FINNED TUBES ----------------- 69 


1. Condensation on Fins in General --------- 69 

2. Significance of Fin Spacing ------------- TS 

3. Significance of Fin Thickness ----------- oo 

4. Significance of Fin Height -------------- Oy 

5. Significance of Tube Diameter ----------- lis 

va CONCLUSIONS AND RECOMMENDATIONS ----------------- ial 

A. CONCLUSIONS 7-9-2997 ee eee oo Ta 

B. RECOMMENDATIONS ------=----=-----------=-----= 120 

APPENDIX A: PREPARATION OF CONDENSER TUBES ----------- NZ 

APPENDIX B: SYSTEM START-UP AND SHUT-DOWN PROCEDURES - 123 
APPENDIX C: RAW DATA OF UNINSTRUMENTED AND 

INSTRUMENTED TUBE ------------------------ 125 

APPENDIX D: RAW DATA OF TUBES TESTED ----------------- Ley 

APPENDIX E: ERROR ANALYSIS -----------------------=---- 150 

LIST OF REFERENCES ------------------------------------ Les 

INITIAL DISTRIBUTION LIST -----93- 9-99-2229 5------------- 6S 


iE. 


fet. 


DW 


Wo. 


ean. 


ele tT . 


eX 


Ze 


a. 


Poll 


ry. 


eo 


evel 


avril. 


LIST OF TABLES 


Geometry of Finned Tubes =----<<---------------- 42 
Instrumented Tube Results for the Vacuum Run --- 61 
Instrumented Tube Results for the 

Eee) hm a a mo 6L 
Constants of Equation (5.6) for All Tested 
a eZ 
Raw Data of Uninstrumented Smooth Tube --------- P25 
Raw Data of Instrumented Tube ------------------ I26 
Raw Data of Tube with Fins of s = 0.5, t = 1.0 

ava ce sre) TET em 128 
Raw Data of Tube with Fins of s= 1.0, t = 1.0 

eee Oe Ze 
Raw Data of Tube with Fins of s = 1.5, t = 1.0 

eT ea a 130 
Raw Data of Tube with Fins of s = 2.0, t = 1.0 

ee ce one 
Raw Data of Tube with Fins of s = 4.0, t = 1.0 

ie a a2 
Raw Data of Tube with Fins of s = 0.5, t = 0.75 

Co. Ot ee eo 
Raw Data of Tube with Fins of s = 1.0, t = 0.75 

and e = 1.0 mm ----999---9 9-H eo ne 134 
Raw Data of Tube with Fins of s = 1.5, t = Q./75 

EN a ESS. 
Raw Data of Tube with Fins of s = 2.0, t = 0.75 

and e = 1.0 mm --------------------------------- 136 
Raw Data of Tube with Fins of s = 4.0, t = 0.75 

and e = 1.0 mm --------------------------------- oe 
Raw Data of Tube with Fins of s = 0.5, 7 = 0.5 

and =e = 170 mm —-====——--——--—-—----—--—-—-~——_—~-~_—__—- Ec 


AE lea 
XIA. 
XX. 
AAT 
otal. 
re oie 
KOE 
DONE 
XAXVI. 


xy, 


POG i 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 
and 


Raw 


Raw 


Data of 
C= c1F0 


Data of 
er= aa 


Data of 
e = 1.0 


Data of 
e-— Jee 


Data of 
e = 2.0 


Data of 
e = 20 


Data of 
ae 0, 


Data of 
e = 2.0 


Data of 
e = 1.0 


Data of 


iL. 2S hotel) e 


Data 


anmeeam @— ep awe a a a ee ee 6 GD eet oe eee eee ee 


rere oe Pe ee se SS owe 8H SPeew si Fs 2 ee SP awe wet SE SP ape Pew SP am am == == 


with Fins 


anmaeear G@aEe @G@eamamwm EG @eamt Gea Ee Ge ae 2 2am 2 4a 42 4a a= 


of Smooth Tube 


fect OF br LGURES 


Cones oeee men elomoes lop lane or Lube —————————— 
Typical Schematic of Retention Angle ------------ 
Schematic of Test Apparatus --------------------- 


Schematic Details of Test Section (Insert 
not Shown) ------------------- ---- -- ---- -- -- - 


Schematic Details of Purge System and Water 
SUMP -- ote oor re rn re 


Photograph of the Instrumented Tube ------------- 


Photograph of a Set of Tested Tubes with 
Fins of t = 0.75 and e = 1.0 mm ----------------- 


Photograph of the Insert Together with a 
Finned Tube Tested ------------------------------ 


Mixing Chamber Calibration Curve ---------------- 


Temperature Distribution for Instrumented 
Tube (Vacuum Run) ------------------------------- 


Temperature Distribution for Instrumented 
Tube (Atmospheric Run) -------------------------- 


Correlation of the Water-Side Coefficient 
Obtained with Instrumented Smooth Tube ---------- 


Comparison of Smooth-Tube Data with Fujil 
Correlation ------------------------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop -------------------------------- 


Variation of Condensing Coefficient with 
Heat Flux ----------------------- ---- ----- ------- 


Cross Plot Showing Best Fin Spacing Among 
Tubes with t = 1.0 and e = 1.0 mm --------------- 


Miniscus Profile for Various Fin Sapcings 
Cyeesce—- leOemm, (b>) s.= 1.5 mm, (c) s = 2.0 mm -- 


2 
24 


yk 


Bez 


oF 


40 


41 


44 


BZ 


2), 


60 


63 


67 


70 


13 


Tad 


719 


he, 


el 


eal Z 


eS 


- 14 


a 


Pn 


sale 


wis 


Sal D 


- 20 


sae 


Normalized Enhancement Ratio Versus Fin 
Spacing for Tubes with Fins of t = 1.0 mm 


and e = 1.0 Wy soe — = == a a ee ee 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 


and e = 1.0 mm (Vac. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for Tubes with t ="1.0 and 


e = 1.0mm (Vac. Runs) ----------~------=-------—- 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 


and e = 1.0 mm (Atm. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


t = 1.0 and e = 1.0 mm (Atm. Runs) ------------- 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.75 


and e = 1.0 mm (Vac. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


t = 0.75 and e = 1.0mm (Vac. Runs) ------------ 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.75 


and e = 1.0 mm (Atm. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


t = 0.75 and e = 1.0 mm (Atm. Runs) ------------ 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.5 


and e = 1.0 mm (Vac. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


t = 0.5 and e- = 1.0 mm (Vac. Runs). —= eee 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.5 


and e = 1.0mm (Atm. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


t = 0.5 and e = 1.0 mm Aen. Runs Se as ae ae 


10 


80 


83 


84 


85 


86 


87 


88 


89 


oO 


oat 


IZ 


99 


94 


22 


720 


24 


as 


26 


we / 


a2 9 


Zo 


- 30 


opi 


Res 2 


3 


34 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 


and e = 2.0 mm (Vac. Runs) --------------------- 


Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 


Pe eee Ct vr mn (Vaen Runs) —--—---------- 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 


and e = 2.0 mm (Atm. Runs) --------------------- 


Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 


and e = 2.0 mm (Atm. Runs) --------------------- 


Cross Plots Showing Best Fin Spacing Among 
Sets of Tubes with Fins of t = 0.5, 0.75 and 


Oe in ieee eee eee a 


Influence of Fin Thickness on Enhancement 
Ratio for Tubes with s = 1.0 mm, for both 


Pressure Conditions --------------------- ------- 


Variation of Condensing Coefficient with 
Heat Flux for Tubes with Fin Geometry Shown 


in the Figure (Vac. Runs) ---------------------- 


Variation of Condensing Coefficient with 
Heat Fillux for Tubes with Fin Geometry Shown 


in the Figure (Atm. Runs) ---------------------- 


Cross Plot Showing Best Fin Spacing for 
Tubes with e = 2.0 mm and t = 1.0 mm, for 


Wah erecoutor Cone Tons. ————————— ———— = —_ _- ___— 


Comparison of Enhancement Ratio Obtained on 
Sets of Tubes with e = 1.0 and e = 2.0 mm 


(Vac. Runs) --------------- -- - --- - - - -- - -- -- 


Condensate Profiles for Tubes with Fins of 


Different Height ------------------------------- 


Normalized Enhancement Ratio of Sets of 


Luce ieeh Sot nO and 2 sO .nm Height ——-—-——-— 


Cross Plots Showing Best Fin Spacing for 


Tubes of Different Diameter -------------------- 


Error Bands for a Tube Tested Under Vacuum 


Gaile takon see ee ee 


Error Bands for a Tube Tested Under 


Atmospheric Conditions ------------------------- 


2) 5) 


iG 


oe 


98 


BK 


eff 


fg 
eee ial 
A 


bist OF sSyiipe ss 


Effective Area of Finned Tube (See kn. 1.10) 
Outside Area of Finned Tube 
Outside Area of Smooth Tube 
Constant Used with Sieder-Tate-type Equation 
Constant Used with Sieder-Tate-type Equation 


Correction Factor for Cooling Water Viscosity 
Changes 


Constant Pressure Specific Heat of Cooling Water 


Equivalent Diameter Used with Beatty and Katz 
Equation 


Diameter to Outside of Fins 
Inside Diameter of Condenser Tube 


Outside Diameter of Smooth Tube, Equal to Root 
Diameter of Finned Tube 


Fin Height 
Z 
eo eee) Gan 
Local Gravitational Acceleration 


Average Heat-transfer Coefficient Given by 
Beatty and Katz Model 


Heat-transfer Coefficient Based on Smooth Tube 
Diameter 


Heat-transfer Coefficient on the Flooded Part 
of the Tube 


Specific Enthalpy of Vaporization 
Water-side Heat-transfer Coefficient 


Effective Heat-transfer Coefficient on the 
Flooded Part of the Tube 


iz 


id) = ake: 
O 


OW 


RW 


Nu 
Nu 


A 


PyarParkP 


eae 


i. 


Re 


Re 


Pap oPop 


Steam-siae Heat-transfer Coefficient 


Average Heat-transfer Coefficient Given by 
Owen et al. Model 


Average Heat-transfer Coefficient Given by Rudy 
and Webb Model 


Heat-transfer Coefficient on the Upnper Tube 
Pert1on 


Thermal Conductivity of Cooling eter 
Effective Thermal Conductivity 

Thermal Conductivity of Condensate 
Thermal Conductivity of Fin Material 
Actual Cooling Water Mass Flow Rate 
Calculated Cooling Water Mass Flow Rate 
Condensing Tube Length 
Logarithmic-mean-temperature Difference 
Fin Density 

Water-side Nusselt Number, nD. /k. 
Steam-side Nusselt Number, NoDo/K¢ 


Pressure Difference Between Points C and B 
Or Poli com wands Ber Gaels. | ) 
Pressure in Condensate (Fig. 1.1) 


Prandtl Number of Cooling Water, C__iu /K . 


pers 
Pressure of Steam 
Heat Flux Based on Outside Area of Smooth Tube 
Coolant-side Reynolds Number, a) aa 
Steam-side Reynolds Number, p-V Di/ues 


ali) (s) orsioha ote) 


Fin Thickness 


iS 


th alae 
Tow 2's 
S 


aenete 


WO 


Cooling Water Inlet Temperature 
Temperature Of Steam 
Cooling water Outlet Temperature 


Average Outer all Temperature, Temperature at 


Root of Fins 
Local Temperature Drop Across Condensate Film 


Average Temperature Drop Across "Gondensa ce saimimm 

(T ae) 
S wO 

Cooling Water Temperature Drop 

Radius of Curvature of the Miniscus Profile in 


the Interfin Space at the Top Part of the Tube 
@nake, “IEF I, 


Overall Heat-transfer Coefficient 
Cooling Water Velocity 

Steam Velocity 
Average Fin Height (Defined by Equation (1.11)) 


Coos natemAsars “(Ff 1p as) 


Coordinate Axis] (mig. ae 


Greek symbols 


Angle Measured from Top of Tube 


Viscosity of Cooling Water at Bulk Temperature 
Viscosity of Condensate at Film Temperature 


Viscosity of Cooling Water at Inner Wall 


Temperature 
Density of Cooling Water 
Density of Condensate 
Density of Steam 


Surface Tension of Condensate 


14 


Piet Cc T1ency 
Condensate Retention Angle Measured from the 


Bottom of the Tube to the Position at Which 
the condensate Just Fills the Interfin Space 


ir 


ACKNOWLEDGEMENTS 


This thesis was completed because of the careful and 


good-willing help of ehe following = peusonce 


sl 


2 


4. 


Dr. P. J. Maneo, siyead sens 

Dr. J. W. Rose, co-advisor in the project. 

Dr. A. S. Wanniarachchi, my second advisor. 

Mr. Tom Christian and the workers in the machine 
shop of the Mechanical Engineering Department, and 


Artemis, my wife. 


I would like to express to them my deep feelings by just 


Saying "Thank you very much." 


16 


Pare nOoOUCT TON 


Modern marine propulsion plants must be designed to 
provide maximum power with small size, low weight and 
reasonable cost. A heavy power plant may have an impact on 
ship dimensions and displacement since the power required to 
achieve a given speed increases with the displacement. A 
sizeable propulsion plant will occupy space that is valua- 
ble, while the machinery room will have poor accessibility. 

A significant decrease in size and weight of a steam 
power plant can be achieved by reducing the size of the 
condenser. Present condensers are large because of the use 
of smooth tubes and generous design margins. The heat- 
transfer performance of any condenser tube is limited by the 
thermal resistances on the water-side, in the wall and on 
the steam-side of the tube. Reduction of any of these 
resistances will contribute to a lighter and smaller 
Gemaenser. 

In recent years, various efforts have been made in order 
to increase the outside (1.e., steam-side) heat-transfer 
coefficient, and these methods have been thoroughly reviewed 
by Marto [Ref. 1]. Outside enhancement techniques include 
Mmicemiuse Of low integral fins, roped tubes, fluted tubes, 
drainage strips attached on the tubes and applied coatings 


to promote dropwise condensation. 


lee 


A. FILMWISE CONDENSATION ON FINNED TUBES 

Filmwise condensation heat transfer on a smooth tube is 
limited by the thickness of the condensate layer. This 
layer forms a relatively high conduction resistance between 
the vapor and the tube surface. To enhance filmwise conden- 
Sation heat transfer, it is then important to thin this 
condensate film. One common way to decrease this condensate 
film is to use radial fins attached to the external surface 
of the tube in order to promote surface-tension effects. 

The theoretical treatment of this problem is very diffi- 
cult owing to the large number of controlling parameters 
Such as the three-dimensional condensate flow pattern, vapor 
shear, surface-tension forces, wall conduction effects and 
condensate retention. Studies following Gregorig [Ref. 2] 
have demonstrated the important role of surface tension in 
modifying the condensate film over the fin surface. The 
presence of fins causes three main effects: 

l. the condensing area increases, 

2. the surface-tension forces thin the condensate film 
over the tip of the fins causing improved heat 
transfer, and 


3. the fins may cause condensate retention especially 
on the lower part of the tube. 


The first two effects are beneficial, while the third 
undermines the heat-transfer performance. The second effect 
is most important in the upper portion of the tube, where 
there is little condensate retention. For steam condensa- 


tion data on finned tubes, Yau et al. [Ref. 3] measured 


3 


Simancement Gatios (1.e€., the ratio o£ vapor-side coeffi- 
cient of a finned tube to the value of smooth tube) in 
excess of the increase in condensing area. This additional 
enhancement 1S most probably due to the thinning effect of 
the surface-tension forces acting on the condensate film. 


This phenomenon can be explained by referring to Figure 1.1. 


Condensate 





Figure 1.1. Condensate Profile at Top Part of Tube 


Since the condensate film hasS a concave shape at point 
A, the pressure within the condensate at this point is 


smaller than the local vapor pressure, and is given by: 


7 pe 
Cen Py = eee) 


> 


A similar expression may be written at point C. 


= eo 
ne z @lks.2) 
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Note that the local radii of curvature at points A and C are 
very large; therefore, the liquid pressure at these points 
nearly equals the vapor pressure. On the other hand, at 
point B the condensate profile has a more concave shape, 
thus resulting in a lower pressure than at either points A 


Oleg: 


- ~ 9 
P,, aie r (Tass 


As a result, there exists a pressure difference from point A 


to point B as given by: 


AP = 5 0. =) ( 1 349) 


AP = o (— = a, (lose 


In general, these pressure differences are positive, causing 
the condensate to be pulled toward point B (which acts as a 
condensate run-off channel), thus thinning the condensate 
film on the fins as well as in the interfin area. This 
thinning effect, therefore, increases the heat-transfer 


performance in excess of the surface area gain. 
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On the other hand, as mentioned above, the third effect 
eee Lins which 1S to retain the condensate, especially on 
the lower part of the tube, creates a relatively thick 
condensate film, causing lower heat transfer. When the fin 
Spacing is reduced, a high surface-tension condensate will 
bridge the fins; therefore, the heat-transfer performance 
may decrease despite that the condensing area has increased. 
Thus, owing to the competing nature of these effects, the 
selection of a best fin geometry needs considerable attention. 

The effect of fin spacing and fluid properties (such as 
Surface tension) on condensate retention has been well 
studied by various investigators [Refs. 4,5,6]. A detailed 
discussion of these studies is presented in the next section. 

One more important effect that 1s common for both smooth 
and finned tubes is the oscillation of the condensate layer 
resulting from the intermittent departure of condensate 
droplets from the lower part of the tube surface. Whena 
droplet departs, the nearby film stretches (becomes thinner) 
and higher heat transfer takes place locally. After the 
droplet has departed, the film grows again and becomes 
thicker. The frequency and amplitude of this oscillation 
are functions of the heat flux, fin geometry and the proper- 
ties of the fluid. To the best of this author's knowledge, 


this transient phenomenon has not yet been investigated. 
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Bie SURVEY OF ANALYTICAL MODELS AND EXPERIMENTAL WORK 

In 1948, Beatty and Katz [Ref. 7] developed a simple 
analytical model to predict the condensing heat-transfer 
coefficient on a horizontal finned tube. Their model was 
Supported by various other researchers, as stated in 
[Ref. 8]. The final form of the Beatty and Katz prediction 


has the following equations: 








3 
ko -(p Oia 1/4 1/4 
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The major fault with their model was the neglect of any 


surface-tension effects. They therefore assumed that the 


Due. 


Mescclenwcaquatlom can be applied to all the external surfaces 
Sele tune (1.e., tip, sides of fins and interfin space). 
According to [Ref. 8 and 9], this model overpredicts the 
condensing coefficient especially for those fluids that have 
mea Surface tensions (for example, water). 

Rudy et al. [Ref. 5] attempted to include the surface- 
tension effect by modifying the original Beatty and Katz 
Petron. FO Obtain this modification, they first investi- 
gated the retention angle problem. They performed static 
(1.e€., with no condensation) and dynamic (1.e., with conden- 
sation) measurements of the retention angle on tubes with 
fin densities from 0.75 to 1.38 fins/mm. They found that 
the retention angle (see Figure 1.2) can be very well pre- 
dicted (within +10 percent) by the following equation: 

a Oomatl - sep Cipede2)) 

£ ig 

The above equation was tested against many fluids such as 
meomrgerant—-il, pentane and water. It is worth mentioning 
Bat Honda et al. [Ref. 6] and Owen et al. [Ref. 8], 
obtained essentially the same equation by using different 
Megeols. Next, Rudy et al. modelled the condensation process 
assuming no heat transfer in the flooded portion of the 


tube. They developed the following prediction: 


Z ee 
Meee letras) Gives) 


ZS, 





of Retention 


Figure 1.2. Yypmealecehemawne 
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where Dex Peerven oy oCuaerTonm(..o9).  Unrtortunately, this 
model tends to underpredict most of the available experimental 
data by as much as about 60 percent [Ref. 9]. 

Owen et al. [Ref. 8] modified the Rudy et al. model by 
including heat transfer from the flooded part of the tube. 
For the fluids that they tested, their model predicted the 
condensing coefficient within +30 percent. The prediction 


that they developed is based on the following equations: 


= a Y 
ho. (1 - 2h +“, (1.14) 
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The above analysis uses Beatty and Katz prediction (equation 
(1.6)) for the unflooded part of the tube. As mentioned 
earlier, the surface-tension effect is not included in equa- 
tion (1.6). However, as discussed in the previous section 
(Section A), the surface-tension forces play an important 
role in filmwise condensation by thinning the condensate 
Film. Therefore, the Owen model needs further improvements 
to account for the surface-tension effect at the unflooded 
part of the tube. 

Recently, Honda et al. [Ref. 9], developed a model to 
predict the condensing coefficient by accounting for surface- 
tension effects and for the non-isothermal behavior of the 
fins. Their model can predict most of the available experi- 
mental data including ll ftluids and 22 tubes ys-hine=2o 
percent. Honda et al. provide a comparison of these four 
models with experimental data. It is worth mentioning that 
the worst disagreement (within +30 percent) between the data 
and their model occurs for the steam data of Rose [Ref. 3]. 
This is probably owing to the high surface tension of water, 
as the surface-tension effects are still inadequately 
accounted for even in this model. Further, it is shown that 
all the models except that of Honda overpredict these steam 
data. In comparison with the other models, the Honda model 
Shows considerable promise. However, it must be noted that 
this model is rather laborious as it requires considerable 


computer time for iterative finite difference techniques. 
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The Honda model may require further improvements to predict 
steam data more accurately. Also, most importantly, a more 
practical model must be developed so that the designers may 
use it with confidence. Such a model may consist of a 
combination of a few correlations. Thus, it is extremely 
important to generate a sufficient amount of experimental 
data covering all important geometric variables (such as fin 
Spacing, fin thickness, fin height, fin shape, tube diameter, 
etc.) for condensation of steam on finned tubes. 

In 1961, Staub [Ref. 10] conducted steam condensation 
experiments with horizontal fluted tubes under vacuum and 
atmospheric conditions. He used a bell-shaped (rather than 
rectangular) fin profile on a root diameter of 14.9 mm. These 
fins had a pitch of 1.0 mm, a thickness of 0.4 mm at the 
root of the fin and a height of 0.5 mm. He reported a vapor- 
side heat-transfer-coefficient enhancement (i.e., the ratio of 
vapor-side coefficient of finned tube to the value of smooth 
tube) of around 3 for both vacuum and atmospheric pressure 
@emaitions. 

At Queen Mary College of London, Yau et al. [Ref. 3] obtained 
steam condensation data at atmospheric pressure on various 
rectangularly-shaped finned tubes having a root diameter of 
ieee. They found an Optimum spacing of 1.5 mm for a fin 
Popermess OL 0.5 mm and a fin height of 1mm. Their optimum 
tube resulted in a vapor-side heat-transfer coefficient 
enhancement of 3 to 4. 

On-going research is underway at the Naval Postgraduate 


School to study the feasibility of using low integral-finned 


2) 


tubes with steam. A test apparatus has been constructcama, 
Krohn [Ref. ll] to systematically study steam condensation 
on a single horizontal tube. The test apparatus was instru- 
mented and tested by Graber [Ref. 12], and meaner system 
improvements for leak tightness were made by Poole [Ref. 13]. 
Poole conducted experiments both under vacuum and at pressures 
near-atmospheric. He used machined, rectangular fins with 

a fin thickness and a fin height of 1 mm on a root diameter 
of 19 mm. He tested a total of six finned tubes with fin 
spacings of 0.5, 1.0, 1.5, 2.0 , 4.0 and 9.0 mm as weliaea. 
smooth tube. His results showed that the best enhancement 
of 3 to 4 occurred for the tube with a fin spacing of 1.5 mm. 
However, Poole had problems with partial dropwise condensa- 
tion since the total filmwise mode lasted only for about 20 
minutes during his experimental runs. Thus, it became 
necessary to repeat his test runs in order to verify his 


preliminary results. 


Cc. ObJECIae is 
The main objectives of this research effort were, there- 
fOne, Le: 
1) establish a preparation procedure that will ensure 
Filmwise condensation on copper finned tubes. 
2) take data on an instrumented smooth tube to determine 
the water-side heat-transfer coefficient with internal 


enhancement. 
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3) take data on a smooth tube to check the data-reduction 
program and to establish a basis for comparison with the 
finned tubes. 

4) take data on tubes with rectangularly shaped fins of 
various fin spacings, fin thicknesses and fin heights to 


obtain a possible optimum fin geometry. 
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II. DESCRIPTION OF APPARATUS 


Ale OY ole OVE Ry ies 

The apparatus used for this research was essentially the 
same used in [Ref. 13], with several small modifications. A 
schematic of the system is shown in Figure 2.1. Steam was 
generated ina 304.8 mm (12 in) diameter Pyrex glass boiler 
using ten 4000 Watt, 480 Volt Watlow immersion heaters. 
Passing through a 304.8 mm (12 in) to 152.4 mm (6 in) 
reducing section, the steam travelled upward through a Pyrex 
glass section 2.44 m (8 ft) in length, around a 180-degree 
bend and back down a straightening section 1.52 m (5 ft) in 
length before entering the stainless-steel test section. 
The condenser tube to be tested was mounted horizontally in 
the test section behind a viewport to permit visual observa- 
tion of the condensing process. Figure 2.2 shows the 
details of the test tube mounted in the test section. 

Steam that did not condense on the test tube passed into 
a stainless steel auxiliary condenser, and all condensate 
was returned via gravity to the boiler. The auxiliary 
condenser was constructed of two 9.5 mm (3/8 in) diameter, 
water-cooled copper lines helically coiled to a height of 
477 Ue le Seeley) ee 

Cooling water for the test tube was provided, on a once- 


through basis, by two centrifugal pumps connected in series. 
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Water to these pumps was provided through a large sump 
capacity about 0.4 fon which was filled with a continuous 
supply of filtered tap water. The water could be throttled 
Mrem Zero to 0.55 l/s (8.8 gpm), giving a maximum water 
velocity of 4.4 m/s (14.4 ft/sec) through the tube. 

A continuous supply of tap water was used for cooling 
the auxiliary condenser. Throttling the flow of tap water 
through the auxiliary condenser was the means used to vary 
the internal pressure of the test apparatus. 

A vacuum pump was used to create vacuum and for removal 
of noncondensing gases from the test section. The vacuum 
pump was in operation throughout the data collection time 


interval. 


beeeeoLolEM INSTRUMENTATION 

The input power to the heaters was changed through a 
panel-mounted potentiometer. A converter with an input 
voltage of 440 VAC generated a signal, which was fed to the 
data acquisition system in order to calculate the power 
input to the heaters. A complete explanation of the con- 
verter function is given by Poole [Ref. 13]. 

The internal pressure of the system was measured manually 
by a U-tube, mercury-in-glass manometer graduated in milli- 
meters. Unavoidably, steam would condense in the manometer. 
Therefore, the varying height of the water column in the 
manometer needed to be accounted for when measuring the 


System pressure. 


So 


All temperatures were measured using 0.25 mm diameter, 
copper-constantan thermocouples. Two thermocouples were used 
for the steam, one for the condensate return and one for the 
ambient. The calibration procedure for these thermocouples 
is described in [{Ref. 13]. The temperature rise of the cooling 
water through the test Section was measured by two Hewleuu. 
Packard (HP) 2804A quartz thermometers as well as by a LOG 
Junction, series-connected, copper—cons tantan tnermepiice 
The temperatures were fed directly into the data-acquisition 
system as described below. During this entire investigation, 
these two independent means of cooling-water-temperature rise 
measurements agreed to within +0.02 K. Note that the coolant 
temperature rise was as low as 0.7 K for the smooth tube 
(under vacuum conditions), while it was in excess of 1720 
for the finned tubes. Thus, the coolant temperature rise was 
measured with an uncertainty better than +3 percent for finned 


tubes. 


C. DATA ACOUISITION TS ore 

An HP 3497A Data Acquisition System, as controlled by 
an HP 9826A computer, was used to monitor the system tempera- 
tures and the converter's signal representing the power input 
to the heaters. The data were immediately processed and also 
were stored on computer disks for later reduction and 


evaluation. 
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Peeeeov oleh MODIFICATIONS 

For some of the early data taken near atmospheric pres- 
Sure, the measured steam-side heat-transfer coefficient, 
when operating with a closed system, disagreed by as much as 
+30 percent between different data runs. This was later 
found to be due to the malfunctioning of the relief valve 
provided on the apparatus for safety reasons (see Figure 
2.1). This valve had a designed relief gauge pressure of 
52 mmHg. However, occasionally it was found to open with a 
system gauge pressure of only about 7 mmHg. Once opened, the 
valve remained open even if the system operating pressure 
was 10 mmHg below atmospheric, thus allowing air to enter 
the apparatus. To avoid this problem, subsequent runs were 
made without the relief valve and a new, more reliable, 
valve was ordered. In the meantime, extreme care had to be 
exercised to avoid system pressurization. | 

Originally, the apparatus was provided with a 10.6 kW 
(3 ton) air-conditioning system to cool the cooling water 
Sump. However, as also mentioned by Poole [Ref. 13], this 
system was inadequate to keep the sump temperature from 
rising. For this reason, this air-conditioning system was 
removed. By removing this unit, the cooling water friction 
losses were reduced and the maximum water velocity was in- 
creased from 4.2 to 4.4 m/s. In the previous stage of the 
investigation, fresh tap water was continuously fed to a sump 


while an equal amount of water was being drained, maintaining 


2 


a constant sump level. Since the sump was small, the hot 
water from the test Section had to be returnedwro the veunee 
For this reason, there was a continuous variation of the 
Sump water temperature. In order to solve this problem, the 
Sump was replaced by an aluminum tank three times its size 
and the outlet hot water was drained directly after its 

exit from the test tube. 

The air ejector used to create vacuum in the test sec- 
tion could not provide a successful operation owing to the 
cycling of the feeding air pressure. Occasionally the pres- 
Sure in the test section could become smaller than the pres- 
Sure in the ejector's air line resulting in a back flow 
which could contaminate the apparatus. This problem was 
solved by the use of a vacuum pump. The suction line for 
this vacuum pump was a penetration in the auxiliary condenser. 
Unavoidably, there was a continuous removal of steam together 
with the non-condensing gases. To avoid steam from reaching 


the vacuum pump (vacuum pumps are designed to remove only 


air), the mixture was passed through a heat exchanger where 
the steam was condensed. This heat exchanger was installed 
as shown in Figure 2.3. Note that no new supply of cooling 


water was used for this heat exchanger; instead, the tap water 

supply line into the water sump was used for this purpose. 
Careful and persistent investigations into the occurrence 

of partial dropwise condensation on the copper tube revealed 


that the presence of nylon components, such as tube bushings 
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(see Figure 2.2) and numerous ferrules used in the apparatus, 
were a source of system contamination. This problem worsened 
when the system pressure was reduced for a vacuum run after 
the system had operated at near-atmospheric pressures. This 
was thought to be owing to the fact that outgaSSing increases 
with increasing temperatures and decreasing pressures. To 
solve these problems, new Teflon tube bushings were manufac- 
tured. Also, all nylon ferrules were replaced with Teflon 
ferrules. Note that Teflon has an outgaSsing rate about one- 


hundred times smaller than nylon [Ref. 14]. 


RE. VACUUM INTEGRITY 

Vacuum tightness for any condensation experiment, 
especially under low pressure conditions, is of vital impor- 
tance. Any data taken while the apparatus is inadequately 
air-tight, will lead to false conclusions because of the 
deleterious effect of the non-condensing gases. Many efforts 
were made in the past to achieve a leak-free system. During 
this thesis effort, there were leak tests on a routine basis. 
A total of eight vacuum-tightness tests were carried out over 
a period of ten months, including the period devoted to taking 
data. The leak rate was always less than that represented 
by a pressure rise equal to 2 mmHg in a 24-hour period. 
These tests were carried out at a pressure of about 85 mmHg 


absolute. 
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feet UB Bo TESTED 

An instrumented tube was fabricated from a thick-walled 
@epeer Cube with an inmaer diameter of 12.7 mm (1/2 in) and 
an outer diameter of 18 mm (3/4 in). The condensing length 
was 13.3 mm. This tube was provided with six wall thermo- 
couples, equally spaced around the tube periphery and 
located radially at the mid-point of the tube wall thickness. 
Details about the fabrication and instrumentation of this 
tube are given by Poole [Ref. 13]. Figure 2.4 shows a photo- 
feopimor the instrumented tube. This tube was used to 
determine the leading constant of a Sieder-Tate type equation 
for the inside heat-transfer coefficient. 

A smooth tube of copper with dimensions similar to the 
Pisce bumented tube was also fabricated. The purpose of the 
smooth tube was to determine the Sieder-Tate leading constant 
Beaeoceeng a modified Wilson plot, and to compare the data for 
the steam-side coefficient with the Nusselt prediction. 

Twenty one finned tubes were also fabricated. These 
tubes had the same inside diameter and length dimensions as 
the smooth tube. Further, they had a root diameter equal to 
the outside diameter of the smooth tube. The fins were 
rectangular in shape. The fin dimensions were systematically 
changed as shown in Table I. Figure 2.5 shows a photograph 
Smee Group of tubes with a fin thickness of 0./5 mm and a 


mie neight of 1.0 mm. 
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2.5. Photograph of a Set of Tested Tubes with 
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TABLE I 


Geometry of Finned Tubes 


Fr aia Pana 

Ban ope cme Thickness Heighw 
(mm) (mm) (mm ) 
0.5, 1.0, 1. Seca 436 La i 
0.5, 1.0, 1. Sie2e oe 0.25 1 26 
0.5, 1:0, 1.5 Zao 40 025 LG 
1.0, 1. 5;02 a0.) 46 an Ze 
1.0 15 1 ae 
1 25 12 Lea 


G. iNSEoRS 

Early investigations showed that with no internal en- 
hancement, the water-side thermal resistance could be as much 
as 90 percent of the overall resistance. Under such condi- 
tions, even a small discrepancy on the water-side resistance 
will translate into a substantially greater discrepancy on 
the condensing heat-transfer coefficient. Therefore, in 
order to improve the accuracy of the condensing heat-transfer 
coefficient, a decision was made to enhance the inside heat- 
transfer coefficient. 

For this purpose, a spiral insert was fabricated. The 
insert consisted of a 6.4 mm diameter stainless-steel rod 
together with a copper wire that was wrapped around the rod. 


The diameter and the pitch of the wrapped wire were 3.2 mm 
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and 20 mm respectively. After wrapping, the wire was properly 
soldered on the rod and then it was machined down to give a 
clearance of 0.5 mm between the outer wire diameter and the 
tube inside wall. This clearance was considered important 

EO avoid heat conduction from the tube wall to the insert. 
Figure 2.6 1s a photograph of the insert with one of the 


finned tubes. 
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III. SYSTEM OPERATION AND PROBLEMS 


Bach tube to be tested was cleaned following the proce- 
dures listed in Appendix A. After installing a tube, the 
system was brought to operating pressure by following the 
procedures listed in Appendix B. The data were collected 
when steady-state conditions were achieved which were 
assumed to occur when the temperature of the steam in the 
test section and the inlet-outlet water temperatures were 
stabilized. Double data sets were taken for each flow 
setting starting with a cooling water flow rate of 80% 

(4.4 m/s). The next water flowmeter settings were 70, 60, 
wae 3, 26, and 20 percent followed by 55 and 80 percent. 
These settings were selected to provide nearly equal-spaced 
heat flux data points. 

Since the tubes were fabricated from copper, a highly 
conductive material, the time interval required to achieve 
steady-state conditions was not more than two minutes. A 
whole run (18 data sets) was completed in a time interval of 
less than one hour. After each change in the cooling water 
flow rate, the system pressure changed slightly. In order 
to keep this system pressure nearly constant, it was neces- 
sary to adjust the flow rate through the auxiliary condenser. 

During the early stages of this experiment, it was 


observed that the thermocouple readings were disturbed by an 
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external radio-frequency (RF) signal. This RF signal was 
generated by the special device used to control the power 
input to the boiler heaters. This device consisted of a 
Silicon-controlled rectifier that chopped the input signal 
to obtain the necessary electrical power, while the chopped 
portion was emitted as an RF signal. This RF signal caused, 
for example, the steam temperature to fluctuate within +4 K, 
while the thermopile readings fluctuated within +0.5 K (note 
that the thermopile readings measuring temperature rise of 
cooling water can be as low as 0.5 K). To alleviate these 
problems, all thermocouples were shielded against the RF 
Signal using braided copper cables. This substantially 
reduced the temperature fluctuations. For example, after 
ene shielding was installed, the thermopile readings fluctu- 
ated less than +0.04 K. Further efforts revealed that the 
small fluctuations which persisted were due to the presence 
of a cable ane carried an indication of the voltage to the 
boiler heaters. This voltage is of "DC" type and it was fed 
to the data acquisition system. Unfortunately, this voltage 
was not truly DC; instead it consisted of the RF disturbances 
described earlier. When this cable was disconnected, abso- 
lutely no fluctuations were observed. Thus, a switch was 
added on this cable to disconnect the source of the RF 
Signal while thermocouple readings were being recorded. Soft- 
ware was modified to gather all thermocouple readings while 


this switch was turned off. When the computer was ready to 
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'read' the voltage to the boiler heaters, the program was 
coded to ‘pause’ and then remind the operator to turn on the 
Switch. Immediately after taking this reading, the computer 
Bemenas Ene Operator to turn ert this switch again. 

The appearance of the condensate film was regularly 
observed through the viewport. In the event the condensate 
film appeared to be patchy, the data run was discontinued 
and the data already collected were discarded. During some 
runs, the film appeared to be totally filmwise, but the 
final data set indicated a condensing coefficient up to 10 
percent greater than the value for the initial data point. 
Recall that the initial and final data sets were for the 
Same cooling water flow rate. It is worth noticing that it 
was easier to obtain filmwise condensation under vacuum 
conditions than at near-atmospheric pressures. This is 
probably due to lower vapor velocities associated with 
atmospheric runs (2 m/s for vacuum and 1 m/s for atmospheric 
conditions). All the data presented in this thesis showed 
less than +5 percent disagreement in condensing coefficient 


Peeveen initial and final data sets. 


A. THE DROPWISE CONDENSATION PROBLEM 

The main objective of this thesis effort was to take 
data under strict filmwise condensation conditions using 
COpper tubes. Pure copper has low wettability with water in 


the presence of a minute contamination; thus, it has the 
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tendency to condense steam under the dropwise mode. Such a 
minute contamination is generally difficult to avoid. A 
slightly hydrophobic sumiace condi tiongcanwheadere. Ver, 
false conclusions since the dropwise condensation mode of 
heat transfer is much more effective than the filmwise mode. 

In the first stage of this thesis effort, the tubes were 
prepared following the procedure recommended by Poole 
(Ref. 13]. This procedure included the cleaning of tubes 
with a mixture of sodium hydroxide and ethyl alcohol (equal 
amounts by mass) at a temperature of about 80 °C. This proce- 
dure waS soon proven inadequate as the copper tube showed 
partial dropwise condensation within 20 minutes. This was 
verified during a four-hour endurance test. The water flow 
rate was held constant and every 30 minutes the temperature 
rise of the cooling water was recorded. This rise continu- 
ously increased and gradually the condensation changed to 
the dropwise mode. 

As a solution to this problem, it was decided to oxidize 
the tube surface more thoroughly using the above-mentioned 
mixture (50 percent) aqueous sodium hydroxide mixed with an 
equal part of ethyl alcohol). The tube was painted with this 
mixture and was then heated with steam for about an hour 
with no water flowing inside the tube. This process resulted 
in the formation of a dark oxide layer on the copper surface, 
which had high wetting characteristics and an insignificant 


thermal resistance. After the oxidation, new endurance tests 
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were made. For a time period of three hours, the cooling 
water temperature rise was essentially unchanged, while the 
condensate film appearance was excellent. Steam-side coeffi- 
cients obtained with and without oxidation agreed to within 
+2 percent. This was probably due to the avoidance of even 
minute areas undergoing dropwise condensation. 

The most significant achievement of this oxidation was 
the reduction of data scattering (less than +l percent), and 


the very good repeatability of similar runs (within +3 percent). 


B. VAPOR VELOCITY LIMITATIONS 

The major assumption in the Nusselt theory is that the 
vapor velocity does not induce shear forces on the vapor- 
liguid interface. This assumption is not realistic in prac- 
tice. In order to regulate the pressure in the test section, 
it is necessary to have some steam condensed in the auxiliary 
condenser. This means that there must be steam flow into 
the auxiliary condenser. The steam flow rate (i.e., the 
vapor velocity) is a function of the power input to the 
heaters. In order to achieve low vapor velocities, the power 
input should be decreased. At low pressures, however, it 
was evident that the fluctuation of system pressure increased 
with a decrease in boiler power. This fluctuation was due 
to the intermittent breakup of vapor bubbles on the water- 
vapor interface in the boiler. Surprisingly, this small 


fluctuation in system pressure, and hence vapor velocity, 
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reflected in the cooling water outlet temperature. To 
minimize this fluctuation, a vapor velocity of approximately 
2m/s past the tube was selected for vacuum runs. For the 
atmospheric runs (maximum power input), the cooling water 
outlet temperature fluctuation was less than +0.005 K, while 
for the vacuum runs (half of maximum power input), the fluc- 
tuation was less than +0.05 K. 

A low vapor velocity is also undesirable from the film- 
wise point of view. When the vapor flow is low, then the 
condensation rate decreases and the film can very easily 
break, resulting in a change of the condensation mode. Once 
the film has broken, there is no possibility for the tube to 
resume the filmwise mode. Also, when the vapor velocity is 
high, there is a very strong wash-off effect on the tube, 
resulting in the reduction in the tendency of the tube to 
change the condensation mode. 

It was decided to use 2 m/s and 1 m/s vapor velocities 
for the vacuum and atmospheric runs respectively. These 
velocities are considered much smaller than the actual ones 
in condensers. However, due to the limited power input, the 
1.0 m/s used for the atmospheric runs was indeed the maximum 
possible. For the vacuum runs, the 2.0 m/s vapor velocity 
was achieved by using half of the available power. Some 
attempts were made to increase the vapor velocity to 4.0 
m/s, but it was considered unsafe due to the violent and 


rapid explosions of the vapor bubbles in the boiler. 
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IV. DATA REDUCTION 


Pimine Enis thesis effort, all the programs, including 
property functions and the calibration curves for the 
cooling water flowmeter and for all thermocouples, were 
carefully checked. Only a very few minor corrections had to 
be made. The programs have the following noteworthy features: 

1. In order to obtain a meaningful mixing-cup temperature 
at the coolant outlet, a mixing chamber was provided at the 
outlet as shown in Figure 2.2. As also discussed by Poole 
[Ref. 13], the considerable pressure drop created by this 
mixing chamber introduced a temperature rise in the cooling 
water. With no steam condensation taking place, the coolant 
temperature rise was measured for various water velocities. 
The results were plotted as shown in Figure 4.1. When taking 
condensation data, the coolant temperature rise was corrected 
by subtracting the temperature rise due to this frictional 
heating at the mixing chamber. For this purpose, the least- 
Squares line shown in Figure 4.1 was used. 

2. The smooth end sections of the tube, which were neces- 
sary to hold the tube in place (see Figure 2.2), were acting 
aS axlal fins. Since copper has a very large thermal 
Sememectivity, axial heat conduction occurred. This heat 
conduction was incorporated into the data reduction process 
by computing the fin efficiency of the smooth end sections 


Oe the tube. 
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3. The tap water temperature changes considerably depend- 
Pmene Seasem: lt iss about £5 °C in winter and about 
22 °C in summer. Because of these changes, it was necessary 
to incorporate water property variations in calibrating the 
flowmeter. For this purpose, the following function was 
used, as done so in [Ref. 15], to calculate a correction 


factor for viscosity changes. 


- = i0sG5=1. 9644 100 \ikees, 25¢lon>) Tan (4.1) 


Mimemvalue of the above function at the time of flowmeter 
mmmpeoratiOon (1.e., at Tin = 22.0 °C) was 0.995. Thus, in 
order to obtain the actual mass flow rate for an inlet 


temperature of Tin, the following equation was used: 


C 
a = 3 E (4.2) 


act “cale 0.995 
The major steps for the data-reduction programs are 


listed below. 


A. PROGRAM SIEDER 
This program calculates the constants of the Sieder-Tate 
type equation from the instrumented-tube data. 
1. Read the cooling water inlet and outlet temperatures 
and compute the average temperature. 
2. Read the six wall thermocouples and compute the 


average wall temperature. 
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3. Convert the water volume flowmeter setting to mass 
flow rate by dividing the volumetric water flow rate with 
the water velocity. 

4. From continuity, compute the cooling water velocity. 
5. Correct the outlet water temperature for the mixing 
chamber effect and then re-compute the average temperature. 
6. Compute the amount of heat transferred to the cooling 

We be le. 

7. Compute the wall thermal resistance and then find the 
inner wall temperature (from step 2, the average temperature 
at the radial location of the thermocouples, is known). 

8. Compute the log-mean-temperature difference. 

9. Compute the efficiency of the left and right tube 
extensions. 

10. Compute the inside heat-transfer coefficient. 

ll. Compute the Nusselt number. 

12. Compute the leading constant of the Sieder-Tate 
coefficient for this Wdatamsec. 

13. Repeat steps 1 to 12 for all the data sets of the run. 

14. Make a least-squares straight line fitting of the 


foOmme 


Ones l/s cl. fea (4.39 
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The C is the unknown leading constant, while B is an addi- 


tional constant found to give an improved fitting. 
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B. MODIFIED WILSON PLOT PROGRAM 
This program calculates the leading constant of the 
Sieder-Tate equation from the smooth-tube data. 

l. Assume a value for the leading constant of the Sieder- 
Tate equation. 

2. Make the steps l, 3, 4 and 5 of previous program 
Pom DER. ” 

3. Compute the heat flux from the specific heat of water, 
the mass flow rate, the cooling water temperature difference 
and the outside tube surface. 

4. Compute the log-mean-temperature difference between 
the steam temperature and the inlet and outlet cooling water 
meMmMocracures. 

5. Compute the overall heat-transfer coefficient. 

6. Assume an outer tube surface Weniser a Cone in order to 
determine the film temperature for the evaluation of conden- 
Sate properties. 

7. Calculate the outside heat-transfer coefficient by 
using the Nusselt equation. 

8. Compute the outer tube surface temperature. 

9. Iterate between steps 7 and 8 until two successive 
outer surface wall temperatures are within +l percent. 

10. Assume a viscosity correction factor for the 
Sieder-Tate equation. 
ll. Compute the fin efficiency of the extending parts of 


the tube. 


Do 


12. Compute the inside heat-transfer coefficient. 

13. Compute the inner surface temperature in order to find 
the film temperature on the coolant side. 

14. Compute the viscosity correction factor and iterate 
between steps 1l through 14 until two successive viscosity 
correction factors are within +l percent. 

15. Compute the new Sieder-Tate coefficient and iterate 
between steps 2 through 15 until two successive Sieder-Tate 
coefficients are within +0.5 percent. 

16. Compute the Sieder-Tate coefficient for all the data 
sets of the run. 

17. Make a least square straight line fitting to establish 


the final leading constant that satisfies all the data sets. 


C. MAIN DATA REDUCTION PROGRAM (DRP4) 
This program was used for computing the condensing heat- 
transfer coefficient for both finned and smooth tubes. 

l. Make steps 1 through 5 of program "SIEDER." 

2. Compute the log-mean temperature difference between 
the steam temperature, and the inlet and outlet cooling 
water temperatures. 

3. Compute the heat flux: 

4, Compute the overall heat-transfer coefficient. 

5. Compute the wall conduction resistance. 

6. Compute the inside heat-transfer coefficient. 

7. Compute the condensing heat-transfer coefficient by 
subtracting the wall and the inside resistances from the 


overall resistance. 
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V. RESULTS AND DISCUSSION 


Numerous data runs were made using he eedueee 
described in Chapter III. Each tube was tested at least 
three times, on different days, in order to obtain repeata- 
ble results. Data were taken for all the tubes at both 
near-atmospheric and vacuum (85 mmHg, 1.7 psia) conditions. 
Complete filmwise condensation was maintained for the data 
runs and the mass concentration of non-condensing gases was 
held between 0.0 and -1.5 percent during all the testings. 
The negative value was indicative of either a slight super- 
heat or an inaccurate manometer reading. It was estimated 
that an error of +0.2 K in steam temperature or an error of 
1.0 mmHg in system pressure translates into approximately 
mz percent Grror in the non-condensing gas concentration. 
Thus, the non-condensing gas concentration was essentially 


zero to within the accuracy of the measurements. 


A. EINSIDE HEAT={TRANSFER COEFFICIENT 
1. Instrumented Tube Results 
The instrumented tube was used to determine the 
inside heat-transfer coefficient mainly owning to the use of 
Miamiisctrt tO Cnhance the heat-transfer performance. By 
the use of the instrumented tube, the wall temperature was 
determined after averaging the six thermocouple readings. 


Therefore, the main assumption in the calculation of the 


ov 


inside heat-transfer coefficient is that the wall tempera- 
ture profile is uniform. However, Figures 5.1 and 5.2 show 
the experimentally-measured temperature profiles around the 
tube, for the maximum and minimum cooling water flow rate 
settings, for the vacuum and near-atmospheric runs, respec- 
tively. It is obvious that the wall temperature profile is 
not uniform; therefore, this direct averaging procedure may 
introduce an uncertainty the magnitude of which is unknown. 
However, Owing to the time constraint, further evaluation 
of this uncertainty was considered to be beyond the scope 
of this investigation. These figures also show least-squares- 
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It is worth noticing that the temperature profile corres- 
ponding to the high heat flux is flatter than the profile 
corresponding to the low heat flux. Also, the temperature 
profile of the near-atmospheric run is flatter than the 
temperature profile of the vacuum run. It is concluded that 
as the heat flux increases, hi increases, hence ho/hi 
decreases and the temperature profile becomes more uniform; 
therefore, the assumption of an isothermal wall becomes more 
reasonable. This is one more reason why the internal 
enhancement (use of insert) improves the accuracy of the 


data. 
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Tables II and III give the values of the amplitude of the 
coSine fit for all data points taken, under both pressure 


©@naltions. 


TABLE II 


Instrumented Tube Results for the Vaccum Run 


Heat Flux ne 

(MW/m2) (kW/m2. K) ho/hi 4 

0.209 10.89 0.301 0.135 
0.207 11.28 0.348 0.142 
0.199 11.45 0.430 0.158 
0), eal 11.61 0.507 0.168 
0.184 11.76 0.567 0.178 
Oelgs 12.05 0.692 0.185 
0.160 12.09 0.832 0.202 

LABOCE L£if 


Instrumented Tube Results for the Atmospheric Run 


Heat Flux ho 


(MW/m2 ) (kW/m? . K) ho/hi a 

0.488 9.18 0.201 0.115 
0.486 9.19 02201 0.119 
0.479 9.25 0.249 0.133 
0.462 9.38 0.289 0.137 
0.447 9.49 0.346 0.149 
0.417 9.75 Oe 4am 0.164 
0.384 10.05 0.585 0.179 


ok 


Also the same tables contain additional information 
about the heat flux, the outside heat-transfer coefficient 
and the ratio of outside to inside heat-transfer coefficients 
fOr all the coolimg water flow rate settings. 

By using the average reading of all thermocouples 
for the wall temperature, the complete form of the Sieder- 
Tate-type equation was determined to be: 

0.8 51/3 (1c) 0.14 


Nu “= 0.0635 Reyer ) + 26.4 (SZ 


ae 


The constant B = 26.4 is an additional constant found to 
improve the fitting. Note that the leading coefficient is 
about 2.5 times greater than the well-known Sieder-Tate 
coefficient of 0.027 [Ref. 16] for plain tubes. This ssi 
Ficant difference is mainly due to the coiled insert used 

to boost the inside heat-transfer coefficient. A discussion 
of the advantages and disadvantages of the coiled insert was 
given in Chapter II, Section G. Figure 5.3 shows the data 
points for both pressure conditions together with a fit of 
equation (5.2) 

Tt is worth mentioning that the use of thermocouples 
to determine the wall temperature may introduce yet another 
uncertainty (in addition to the incorrect averaging mentioned 
earlier) in the measurement of each local wall temperature. 
This is due to the completely different thermal conductivity 


of the wall area in the immediate vicinity of the thermocouple 
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hot junction. For example, the thermocouple wires are 
covered with plastic; therefore, they present a low thermal 
conductivity. Also, in order to embed the thermocouple, 
there must be sufficient clearance, which certainly presents 
an additional thermal resistance to the heat flow. 

The advantages of the use of the instrumented tube 
are: 

l. It does not depend on other theoretical formulations 
that may have unavoidable errors owing to the assumptions 
made to formulate the problem (for example, the Modified 
Wilson Plot technique uses the Nusselt-type equation for the 
outside heat-transfer coefficient). 

2. It is not too sensitive to vapor shear effects since 
these effects are lumped into the thermocouple readings 
(1.e., the vapor shear reduces the thickness of the conden- 
sate layer; therefore, the thermocouples would give higher 
temperatures). 

A computer program, named "SIEDER," was written to 
reduce the data taken on the instrumented tube. For com- 
pleteness, the raw data taken on the instrumented tube for 
both pressure conditions are given in Appendix C, Table VI. 
Hereinafter, the data reduction method that uses the instru- 
mented tube water-side heat-transfer coefficient will be 
referred to as the "Direct" method. 

2. Uninstrumented Tube Results 
The uninstrumented smooth-tube data were used to 


determine the inside heat-transfer coefficient by the use of 
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the modified Wilson plot technigue. Hereinafter, this tech- 
nique will be referred to as the "Wilson" method. Also the 
smooth tube results concerning the outside heat-transfer 
coefficient were used as a basis for comparison with the 
finned-tube results. 

The analysis procedures used for the Wilson method 
is the same as that given in [Ref. 17]. While a Sieder- 
Tate-type equation was used to represent the inside heat- 
transfer coefficient, a Nusselt type expression was used for 
the Outside coefficient as given by: 

k2o-(p--p.)h, 1/3 
hy = a(==—5- 59) (5.3) 
Hee) 

where a 1s a disposable constant that must be found itera- 
tively. Note that with zero vapor velocity, the Nusselt 
prediction results ina value of 0.655 for a. However, 
Since it is impossible to obtain zero vapor velocity, a 
higher value is generally obtained. During this work, an 
approximate value of 0.72 was obtained for a. 

A computer program, named "WILSON," was written to 
run the Wilson method on the smooth-tube data. A total of 
four runs (two each under vacuum and at near-atmospheric 
pressures) were made. These runs resulted in values of the 
Sieder-Tate coefficient of 0.0716, 0.0714, 0.0700 and 
0.0708, with an average value of 0.071. This value is about 


10 percent higher than the value found with the instrumented 
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tube; 1i.e., the Wilson method results ina 10 percent 
smaller water-side resistance (less conservative) than the 
resistance predicted by the use of the instrumented tube. 
Again for completeness, the raw data taken on the uninstru- 
mented tube under both pressure conditions are given in 
Appendix C, Table V. 

3. .Summary Of the Inside sheat=lmansicms eee em menage 

In order to compare the two above-mentioned methods 

(Direct and Wilson), the data taken on a smooth tube under 
vacuum conditions were reduced to yield the outside heat- 
transfer coefficients using the water-side coefficients 
derived by each method. For this purpose, these data were 
plotted in Figure 5.4 together with the original Nusselt 
eguation for a single horizontal tube, re-written in terms 
of the indicated non-dimensionalized terms: 


Nu 
S _ io 
=. 5 Ome 5 ar 


Re 


(5.4) 


Also shown in this figure is the correlation of Fujii, et 


al [Ref. 18]: 


—-S_ = 9,96 Fl’? 


The Fujii et al. correlation is based on numerous data 


taken on tube bundles in both in-line and staggered arrange- 


ments. It can be seen that the coefficient on the right-hand 
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side of equation (5.5) is larger than the value in equation 
(5.4) owing to the presence of vapor shear (recall that the 
Nusselt equation is derived for zero-vapor-shear conditions). 
For completeness, the error bands of the data used for this 
correlation are also shown. 

Note in Figure 5.4 that the data of this investiga- 
tion on the smooth tube are slightly higher than those 
predicted by equation (5.5). Also, it is worth noting aie 
the use of the Wilson method results in data that are in 
better agreement with the Fujii correlation than the data 
given by the use of the Direct method. Thus, with an appro- 
priate function for the outside heat-transfer coefficient 
and accurate smooth-tube data, the modified Wilson plot 
technique appears to offer a satisfactory and simple method 
for obtaining the inside heat-transfer correlation. However, 
Since the error band of the Fujii correlation is much wider 
than the disagreement between this correlation and the 
present data, it may be unwise to draw any conclusions as 
to the accuracy of these two methods (Direct or Wilson). 

As a result, it is not possible to say which method of measur- 
ing the inside heat-transfer coefficient is more accurate. 
However, except where noted below, all the outside heat- 
transfer coefficients obtained in this thesis were calcu- 


lated using the Direct method. 
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B. PERFORMANCE OF FINNED TUBES 

Data were taken for the finned tubes given by Table I, 
under vacuum (85 mmHg) and at near-atmospheric pressures 
with vapor velocities of 2 m/s and 1 m/s respectively. The 
computed steam-side coefficient was based on the smooth tube 
outside area (1.e., the area if the fins were machined off). 

The data reductions were performed using the program 
named "DRP4." This is a revised version of the program 
(DRP2) used by Poole [Ref. 13]. The raw data taken on finned 
tubes are presented in Appendix D. 

1. Condensation on Fins in General 

Figure 5.5 shows the variation of heat flux with the 

steam-side temperature drop for a set of five finned tubes 
under vacuum conditions. Both the thickness and height of 
the fins were 1.0 mm, while the spacings were 0.5, 1.0, 1.5, 
2.0 and 4.0 mm. Notice that the slope of the curves in- 
creases with fin spacing up to the point where the optimum 
Spacing (1.5 mm) is reached and then the slope decreases 
with further increase in fin spacing. This means that the 
enhancement of the condensation process is highly dependent 
on the fin spacing. The data were fitted with least-squares 


curves using the following relation: 


Cl cla ee (S705 


and AT was computed by using the following relation: 
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g = h AT Cos 


mMicGgend was measured experimentally, while ho was obtained 
by subtracting inside and wall resistances from the measured 
overall resistance. From eguation (5.6), it is possible to 
Calculate the heat flux for any temperature difference be- 
tween the steam and the tube wall when the constants a and 

b are known. 

Table IV gives the values of a and b for all the 
tubes tested for both pressure conditions. Notice that the 
value of the constant a has almost the same range for all 
the finned tested tubes for both pressure conditions. The 
values of both constants a and b increase as the fin spacing 
increases up to the spacing of 1.5 mm (it will be shown 
later that this is the optimum spacing) and then they 
decrease. The range of the exponent b is different between 
the vacuum and near-atmospheric runs; i1.e., the range for 
the vacuum runs is 0.76 to 0.89, while for the near-atmospheric 
mst 15 0.87 to 0.94. 

Figure 5.6 shows a relative plot of the outside 
condensing coefficient for the same set of tubes just 
described. For comparison purposes, the smooth-tube data 
and the Nusselt prediction are plotted on the same figure. 
All tubes tested show good enhancement over the smooth tube. 
It is worth noticing that even the smooth tube results are 


higher than the Nusselt prediction; this is mainly 
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attributable to Nusselt's conservative assumption which 
neglects vapor shear. 
Combining equations (5.6) and (5.7), nye can be 


written as: 


(5 ay 


It appears that equation (5.8) results in curves almost 
parallel to the Nusselt prediction, indicating that the 
effect of heat flux is the same on the finned surface as on 
the smooth tube. 

On the above figure, the dotted lines represent the 
fittings of the data when the data were processed by the 
Wilson method. Notice that the condensing coefficients 
given by the Wilson method are lower by about 10 percent 
than those given by the Direct method (recall that the 
water-Side resistance given by the Direct method is higher 
than that given by the Wilson method). 

It is also shown in the above figure that as the 
heat flux decreases, the vapor-side heat-transfer coeffi- 
cient increases. Note that there is about a 10 percent 
difference between the high and the low heat flux results. 
The reason for this increase is that, as the flux decreases, 
the condensate thickness decreases, resulting in less 


conduction resistance and improved heat-transfer coefficient. 
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Pea omni tl cance Of Fim Spacing 

This section presents results showing the variation 
of the condensing coefficient with heat flux, as well as the 
variation of heat flux with steam-Side temperature drop for 
both pressure conditions, for the four sets of finned tubes. 
Figures 5.10 to 5.25 show that for all the fin thicknesses 
(i.e., 1.0, 0.75 and 0.5 mm) and for each of the fin heights 
(i.e., 1.0 and 2.0 mm), the best performance was obtained 
with the tubes of 1.5 mm fin Spacing. This was true for 
both pressure conditions. The second best performance was 
given by the tubes with spacings of either 1.0 mm or 2.0 mn. 
The third best performance was given by the tubes with fin 
Spacing of 4.0 mm and the tubes with a fin spacing of 0.5 mm 
showed the poorest performance. The performances of the 
last two tubes however were interchanged as the pressure 
conditions were changed from vacuum to near-atmospheric. An 
effort to explain this behavior will be made in the next 
sections, where the influence of the fin thickness and 
height on the heat-transfer performance is discussed. 

As shown in Appendix E (Error Analysis), the error 
in the calculation of the condensing coefficient increases 
as the heat flux decreases. For this reason, the comparison 
of finned tubes should be done in the regime where the error 
is the minimum possible. This means that the data points 
corresponding to the high-heat-flux settings should be 


chosen. The comparison of the finned tubes is made through 
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the enhancement ratio, which is defined to be the tae onem 
the outside heat-transfer coefficient of a finned tube to 
the smooth-tube value (both evaluated on the same heat 
flux). However, the capability of the smooth tube to remove 
heat is limited (see Figure 5.6 showing that the maximum 
heat flux on the smooth tube is only about 30 percent of the 
maximum heat flux on the best finned tube). Thus, it was 
necessary to extrapolate the smooth-tube data to higher heat 
fluxes in order to make the comparison of the finned tubes 
on an equal heat flux basis. 

Choosing 0.35 and 1.0 MW /m as reference heat fluxes 
for the vacuum and near-atmospnerie runs, respecery— a 
cross plots of the enhancement ratio versus the fin spacing 
are given in Figure 5.7 for the set of tubes with 1.0 mm fin 
thickness and height. It is shown that the optimum fin 
Spacing found is 1.5 mm. The 0.5-mm fin spacing which pro= 
vides the most area increment has one of the worst performances. 
This poor performance can be easily explained by the fact 
that the entire interfin area was covered by the condensate. 
Visual observations showed that the condensate retention 
angle was 180 degrees, which is also very well predicted by 
equation (1.12). The thermal resistance induced by this 
thick condensate layer deteriorates any benefit from in- 
creased area; therefore, the heat-transfer performance is 


highly reduced. 
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As the fin spacing increased from 0.5 to 1.0 mn, 
the retention angle decreased and a miniscus in the interfin 
area at the top of the tube appeared (see Figure 5.8(a)). 

For the 1.0 mm fin spacing, the retention angle was about 

110 degrees. The change of the flooding characteristics is 
the reason for the sudden increase of the enhancement ratio 
as the fin spacing changed from 0.5 to 1 mm. Visual obser- 
vations showed that the optimum fin spacing (1.5 mm) resulted 
in a condensate film profile (in the interfin area at the 
top part of the tube) similar to that shown in Figure S36 
Notice that the bottom of the miniscus is just touching the 
tube surface. When the fin spacing is further increased, 

the miniscus radius becomes larger and the retention angle 
reduces even more. At the 2.0 mm fin spacing, the miniscus 
is no longer an arc of a circle (see Figure 5.8(c)). 

As mentioned in the beginning of this chapter, the 
condensing-coefficient calculation is based on the smooth 
tube external diameter; therefore, the effect of area addi- 
tion is lumped together with all the other factors that may 
influence the heat-transfer process. In order to give some 
insight to the effect of all the other factors, the enhance- 
ment ratio was divided by the area gain (i.e., the ratio of 
the condensing area of finned tube to the area of smooth 
tube) of each finned tube in order to obtain the normalized 
enhancement ratio. Figure 5.9 shows this normalized enhance- 


ment ratio. It is clear that the enhancement is due to 
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more than just an area gain. This verifies the measure- 
ments of Yau et al. [Ref. 3], and an explanation for this 
behavior was presented in Chapter I. Also, Figure 5.9 shows 
that the normalized enhancement ratio increases as the 
Spacing increases and that the rate of this effect gets 
smaller for fin spacings larger than 1.5 mm. This reduction 
is attributable to the fact that the pulling action (thinning 
effect) diminishes as the fins get further apart. Equation 
(5.9) gives that the pressure gradient between points C and 

B (see Figure 1.1) is inversely proportional to the increase 


Sein spacing (s). 





ax | - S (5.9) 
CB 


Therefore, the thinning effect decreases as the fin spacing 


increases. It can be concluded that any assumption that 
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Careful observation of Figure 5.7 shows that the 
Per@mermance of the tubes of 0.5 and 4.0-mm fin spacing 
interchanges between the vacuum and atmospheric runs, 
respectively. The heat that was removed during the vacuum 
runs was about three times smaller than the heat removed 
during the near-atmospheric runs; thus, the rate of conden- 
Sate generation for the vacuum runs was three times smaller 


than for atmospheric runs. For all the vacuum runs, no 
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matter which tube was used, the condensate flow pattern was 
invisible. During the near-atmospheric runs, the flow 
pattern was also invisible except for the tubes of 0.5 mm 
fin spacing. In these cases (0.5 mm fin spacing), the 
draining rate was very fast and ripples in the interfin 
Space were clearly visible. If this is the case, then the 
convective action of the condensate layer should enhance 
the condensation process. After this observation, it became 
clear why the performance of the tubes with 0.5 mm spacing 
became better than the performance of the tubes with the 
4.0 mm spacing, for the near-atmospheric runs. 

As mentioned earlier, Honda and Nozu [Ref. 9] 
present a detailed comparison of various theoretical models 
to predict the condensing coefficient on finned tubes. All 
the models have the assumption that the inertia pecan the 
momentum equation and the convective term in the energy 
equation can be neglected in comparison with the other 
terms. The author feels that the above assumption is 
realistic for all the fin spacings except those that corres- 
pond to completely flooded tubes. To further support this 
argument, the reader is referred to [Ref. 9] where comparison 
of the four models is presented with various experimental 
data, among which are the steam data given by Rose [Ref. 3]. 
Honda and Nozu show that the maximum disagreement between 
the data of Rose and predictions (all four models) is for 


the fin spacing of 0.5 mm (completely flooded tube). 
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SB. Sein Leanoe of Fin Thickness 


Three sets of finned tubes were manufactured with 
mene thicknesses of 0.5, 0.75 and 1.0 mm, while the fin 
height was held constant at 1.0 mm. Each set consisted of 
Beivye tubes with fin spacings of 0.5, 1.0, 1.5, 2.0 and 4.0 
mm. The relative plots of the outside heat-transfer coeffi- 
cient for each group of tubes, toate under vacuum and 
atmospheric conditions, were already shown in Figures 5.10 
to 5.21. Notice that the shape of the curves to fit the 
data points is very nearly independent of the fin thickness. 
Figure 5.26 is a cross plot showing the enhancement ratio as 
a function of the fin spacing for both pressure conditions, 
with the fin thickness as a parameter. Notice that for 
almost all the fin spacings, the 0.5-mm fin thickness gave 
the smallest gain even though it provides the maximum 
condensing area. An effort to explain this phenomenon will 
be made at the end of this section. Between the 0.75 and 
1.0-mm fin thicknesses, it is not very clear which of the 
[myonis Superior. For the vacuum runs with fin spacing 
around the optimum value (i.e., 0.5 mm < s < 1.5 mm) it 
appears that the 1.0-mm thickness is a stronger candidate 
for the optimum fin thickness. On the other hand, for 
atmospheric runs with most fin spacings, the 0.75-mm fin 
thickness appears to be the best candidate. 

In order to further study the dependence of fin 


thickness on the enhancement ratio, data shown in Figure 
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5.26 were cross plotted as shown in Figure 5.27. As mentioned 
earlier, initially there were only three tubes with fin 
mmmekmess Of 0.5, 0.75 and 1.0 mm. For vacuum conditions, 
the inadequacy of this fin thickness range can be seen 

since the optimum thickness has not been bracketed. There- 
fore, an additional tube with a fin thickness of 1.5 mn, 

a fin height of 1.0 mm and a fin spacing of 1.0 mm, was 
manufactured. The enhancement ratio found on this tube is 
also plotted in Figure 5.27. It can be seen that the tube 
with a 1.0-mm fin thickness is still the best candidate for 
the vacuum runs. This figure also shows possible error 
bands for the computed enhancement ratios. Since the 
possible error bands are in excess of the enhancement ratios 
at these two fin thicknesses, no attempts were made to 
explain why the optimum fin thickness is affected by the 
pressure condition. Therefore, it can be concluded that the 
Seermum £inm thickness is around 0.75 to 1.0 mm for both 
pressure conditions. 

As the fin thickness increases with fin spacing held 
constant, Ae fin density (fins/mm) decreases. In the 
limit, the tube will become a smooth tube with an outer 
diameter gual to the fin diameter. Recall that in the 
Nusselt equation (5.3), the tube diameter is in the denominator; 
therefore, the condensing coefficient is inversely propor- 


1/4 


tional to the D. for constant wall temperature, or, at 


a constant wall heat flux, the condensing coefficient is 


eS 


inversely proportional to De Thus, 
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where es is based on the outer diameter of the tube, 
including the fin height. Now, if this condensing coeffi- 
cient is based on root diameter (as it is done throughout 


mms thesis), then 


D = fh 


Bose... De (5.11) 


Miemecoloimatlon Of equations (5.10) and (5.11) results in: 


273 
(Seek 2 ) 


This means that there will be very little enhancement by 


increasing the fin thickness to the limit. For example, 


is 
wemlad equal 1.0/7. 


Tee. = 21.0 mm and D = 19.0 mm, the ratio h via 
O EO Ae 


When the fin thickness 1s reduced, there is a finite 
meen Of Lins that can be added in a given tube length 
(i.e€., it is not physically possible to take a smooth tube 
when fins of infinitesimal thickness are present). Since 
there is a Significant increase in heat-transfer area as the 
fins become thinner, an enhancement for the condensing 
process may be expected on the basis of area increment. 
Unfortunately, this was not the case as discussed below. 


Equation (1.12) predicts the retention angle and it is shown 


FOS 


that the condensate retention is independent of fin thick- 
ness. This was visually verified during this investigation. 
Comparing two tubes which have the same length and fin 
spacing, the one that has thinner fins would hold more 
condensate when the retention angle is the same. This means 
that part of the metal of a thick fin is replaced by conden- 
sate, or in other words, the thin-finned tube presents a 
smaller effective thermal conductivity than the thick-finned 
tube resulting in lower performance. Also, as the fin 
thickness decreases, the fin efficiency may play a role for 
tubes made of low-thermal-conductivity metal. However, for 
copper tubes with fins of 0.5 mm thickness and 1.0 mm 
height, the fin efficiency was estimated as 99 percene 
when the outside heat-transfer coefficient was 40 Kh moe 
Since both fin thickness extremes show that the 
heat-transfer process is not enhanced, it appears that the 
optimum thickness found (around 0.75 to 1.0 mm) is reasonable. 
Careful observation of Figure 5.26 reveals that the 
tubes with 0.5-mm fin spacing do not follow the trend of the 
tubes with larger fin spacings. In this case, the maximum 
performance was given by the tube with the 0.75-mm thick- 
ness, then by the tube with 0.5-mm thickness and finally by 
the tube with 1.0-mm fin thickness. This sequence holds for 
both pressure conditions. The rate of condensate flow | 
around the tube was observed to be dependent on the fin 


thickness. The fastest flow pattern was observed for the 
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tube with 0.5 mm fin thickness, then for the tube with 0.75 
mm and finally for the tube with 1.0 mm fin thickness. This 
Surprising observation can be explained if the fin tip 
effect is introduced; i.e., since the fin tip is rectangu- 
lar, there is a possibility that some of the condensate 

will fall into the interfin spacing and some will follow the 
circumferential path. If this condensate flow is realistic, 
it can be concluded that the thicker the fin, the more 
condensate will flow circumferentially along the fin tip. 

To further support this conclusion, numerous runs 
were made at a different fin spacing (1.0 mm) where the 
miniscus in the interfin area at the top part of the tube 
1s visible (recall that the tubes with 0.5-mm fin spacing 
are completely flooded). If the miniscus radius is different 
as the fin thickness changes, it can be concluded that the 
condensate flow pattern is influenced by the fin thickness. 
The tube with 0.5 mm fin thickness had a smaller miniscus 
radius than the tube with 1.0 and 1.5 mm fin thickness. 

This means that the previous conclusion may be correct; 
therefore, when the fin tip is wide and of course rectangu- 
lar, there will be circumferential drainage along the fin 
tip. On the other hand, for tubes with narrow fin tips, 
most of the condensate will fall into the interfin spacing, 
resulting in a very fast flow pattern and the creation of 
ripples. If this is the case, the thin fins should show 


an increased performance due to the convective action of 


Or 


the condensate. However, the final performance is limited 
by the excessive condensate amount that is attached continu- 
ously on the tube. These may be the reasons why the tube 
with 0.5 mm fin thickness has better performance than the 
tube with 1.0 mm fin thickness when the fin spacing is 
0.5 mm. Between the tubes with 0.5 and 0.75 mm fin thick- 
nesses, there is no clear visible difference in the flow 
pattern; therefore, the last tube is superior because it 
presents less conduction resistance (it holds less conden- 
sate), and it certainly has a higher fin efficiency. 

In general, it can be concluded that for fins of 1.0 


mm height, the enhancement ratio was weakly dependent On @iam 
thickness. The most significant reason is that the fin 
thickness has no influence on the retention problem as the 
fin spacing does. 

In order to check the validity of some of the 
previous conclusions, data were taken for one more tube. 
This tube had a fin spacing, thickness and height of 1.25 
mm, 1.25 mm and 1 mm, respectively, and it was made from the 
tube with 1.5 mm fin thickness whose fins were machined down 
by 0.25 mm. A tube with 1.25-mm fin spacing was expected 
to have a better heat-transfer performance than the tube of 
1.0 mm fin Spacing since the condensing heat-transfer coeffi- 
cient increases up to a fin spacing of around 1.5 mn. 


However, the performance should be limited owing to the 


fact that the fin thickness is outside the optimum range 


KG 


0m >) coe. O Tm). Figures 5.28 and 5.29 show the relative 
mers Greene Condensing CcoefEIcient with heat flux for both 
pressure conditions, for three tubes whose fin geometry is 
shown in these figures. It is shown that the performance 
of the new tube is almost the same as with the tube with 
1.0 mm fin spacing and thickness. This means that what was 
gained by the spacing increment (from 1.0 to 1.25 mm) was 
Meee by the fin thickness increase (from 1.0 to 1.25 mm); 
therefore the previous conclusions about the fin spacing and 
thickness effects seem to be reasonable. 
Peon etcanecgon bin Height 

As shown in Table I, two fin heights (1.0 and 2.0 
mm) were examined during this thesis effort. The geometry 
of tubes tested is given in Table I. The heat-transfer 
performance of these tubes has already been shown with 
Preures 5.22 to 5.25, and for both sets of tubes, the 
optimum fin spacing was determined to be 1.5 mm. Figure 
5.30 shows the enhancement ratio for the set of tubes with 
2.Q0-mm fin height. It is shown that the maximum enhancement 
is about 5 and 6.5 for the vacuum and atmospheric runs, 
respectively. Figure 5.3l is a cross plot showing the 
enhancement ratio versus the fin spacing for both pressure 
conditions, with the fin height as a parameter. It is shown 
that the set of tubes with 2 mm fin height has better 
performance than the set of tubes with 1.0 mm fin height. 


The reasoning is two fold: 
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1. Surface area increases with increasing fin height 
(see Figure 5.30). 

2. The retention angle decreases and hence heat-transfer 
performance increases, with increasing fin height as pre- 
dicted by eguation (1.12). This was also visually veri- 
fied during this thesis effort. 

Examining carefully the above mentioned Figure 5.31, 
it can be seen that the performance sequence is not the same 
for both sets of tubes. The tube with 1.0-mm fin Spacing 
is the second in performance among the set of tubes with 
1.Q0-mm height. However, for the set of tubes with 2.0 im 
height, the tube with 1.0-mm fin spacing is the Ehiran am 
performance while the tube of 2.0-mm fin spacing becomes 
second. This observation holds for both pressure conditions. 
This means that as the fin height increases, there is a 
tendency for the optimum spacing to moveto a higher value. 
This tendency is mainly due to the increased amount of | 
condensate which is retained between the fins at the top 
part of the tube, caused by the increased condensate 
generation. This fact was verified visually during the 
experiments, and Figure 5.32 shows a schematic of the con- 
densate appearance on those two tubes which differ only 
by the fin height. It is shown that the interfin space at 
the top part of the tube with 1.0-mm fin height ‘holds' less 
condensate than the tube with 2.0-mm fin height. This obser- 


vation may reveal that the mechanism of thinning the film 
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diminishes when the fin height increases. The pressure 
gradient in the condensate film between points A and B 


(see Figure 1.1) can be approximated as: 





(5. dee 


Thus, the pressure gradient between points A and B decreases 
with increasing fin height (e); therefore, it should be 
expected that the tubes with 2.0-mm fin height would provide 
a less effective thinning effect. To prove this statement, 
one more cross plot was prepared by plotting the normalized 
enhancement ratio versus the fin spacing (1.e., separate 

the area addition effect). Figure 5.33 shows clearly that 
the set of tubes with 2.0-mm fin height is not improving the 
condensation enhancement mechanisms as effectively as ‘the 
set of tubes with 1.0-mm fin height. 

It 1s worth mentioning that as the fin height 
increases, the portion of the fin that extends above the top 
flooded part of the tube increases (see Figure 5.32 and 
compare the extending fin part for the tubes with 1.0 and 
2.0 mm fin height); therefore, this part of the fin presents 
the highest heat-transfer performance since the condensate 
layer around the fin is extremely small. Of course, the 
performance of the tip is limited by the fin efficiency 
which decreases with increasing fin height; therefore, the 
set of tubes with 2.0-mm fin height has better performance 


mainly because of the excessive area addition. 
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Up to now, two ways to increase the external surface 
area of a tube were mentioned; 1.e., either by reducing the 
fin thickness or by increasing the fin height. In the 
previous section it was shown that by changing the fin 
thickness there was not much enhancement, the reason being 
that the fin thickness does not improve the flooding problem 
and that the area addition happens in the regimes of the 


tube that are flooded. On the other hand, by increasingmeee 


0 mm a aim ee 


SSS nae 


transfer performance. 
J. Cle ni PicGarece vor ecm ec emer 

Yau et al. [Ref. 3] took data for filmwise conden- 
Sation on copper finned tubes under near-atmospheric pres- 
sure. The fin shape was rectangular and the fins were 0.5 
mm in thickness, 1.0 mm in height, and the fin spacings were 
from 0.5 up to 20.5 mm (13 different fin spacings). The 
diameter of the tubes (at the root of fins) was 12.7 mm and 
the vapor velocity was 1.1 m/s. 

Figure 5.34 shows the enhancement ratio versus fin 
Spacing for the data taken by Yau et al. and the data taken 
during this thesis effort. The only difference between the 
two sets of tubes was the root diameter, which was 12.7 and 
19 mm for Yau et al. and the present work, respectively. 
It can be seen that there is a similarity in the enhancement 


ratio curves and for both sets of tubes the best fin spacing 
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found was 1.5 mm. However, there iS a significant differ- 
ence in the enhancement ratio, which can be explained by 
the following reasons: 

1. The retention angle increases as the tube diameter 
decreases, as it is predicted by eguation (1.12), resulting 
in a lower condensing coefficient. 

2. The smooth-tube performance is predicted by the 
Nusselt equation, and it gives that a tube with smaller 
diameter will have higher condensing coefficient (recall 
that the tube diameter is in the denominator). 

Therefore, when the enhancement ratio is computed, 
the tube with smaller diameter will have the numerator small 
and the denominator large, resulting in a smaller ratio 


than that which was obtained during this thesis effort. 


ie 


Vic ONCiUS TONS Shoe e Cc OMMENDATIONS 


fee CONCLUSIONS 

1. For the best fin spacing of 1.5 mm, a fin thickness 
of 1.0 mm and a fin height of 2.0 mm, heat-transfer enhance- 
ments around 4 to 5 were realized for the vacuum runs and 
around 6 for the atmospheric runs, while the area increment 
was only about 3 times the smooth-tube area. The additional 
enhancement is due to the thinning effect; therefore, any 
theoretical model that neglects this effect may have severe 
deviations from the experimental data. 

2. The 1.5 mm fin spacing appears to be the optimum 
among the fin spacings tested (0.5, 1.0, 1.5, 2.0 and 4.0 mm) 
for both vacuum and near-atmospheric pressure conditions. 

3. The condensation enhancement does not strongly depend 
on the fin thickness. However, the tubes with fin thickness 
around 0.75 to 1.0 mm performed better than the tubes with 
0.5 and 1.5 mm fin thicknesses. 

4. As the fin height increased from 1.0 to 2.0 mn, 
the condensing coefficient increased only by about 20 percent, 
despite an area increase of about 50 percent. This lower- 
than-expected performance is believed to be owing to the 
poorer thinning effect of the condensate film with increas- 
ing fin height. 

5. As the tube diameter increases, the retention angle 


decreases resulting ina better heat-transfer performance; 


le 


therefore, the tubes with a root diameter of 19.0 mm showed 
a superior enhancement over the tubes with a 12.7 mm root 
diameter. 

6. The consistency of equation (1.12) that predicts the 
retention angle was visually verified throughout all the 
experiments. Certainly, it must be a very useful equation 
to predict the condensate retention angle when modelling 
Steam condensation on finned tubes. 

7. The leading coefficient of the Sieder-Tate-type 
eguation used:to determine the water-side resistance with 
the instrumented tube was 0.0635, while the modified Wilson 
plot technique gave the same coefficient as 0.071. 

8. The test apparatus successfully operates under any 
pressure condition. An extremely tight apparatus is 
achieved resulting in minimization of the non-condensing 
gas effect. 

SinerA tube-cleaning and preparation procedure has been 


established to ensure filmwise condensation on copper tubes. 


B,  eRECOMMENDA EOS 

1. Calibrate the existing pressure transducer in order 
to replace the U-tube, mercury-in-glass manometer. This 
will make each run faster and more automatic (less data to 
be entered manually). 

2. Install the relief valve on the test condenser. 

3. Take data with different shaped fins (at the tip as 
well as at the root of the fin) to investigate the effect 


of fin geometry. 


JL AO) 


4. Take data with tubes that are fabricated with 
different thermal conductivity materials to investigate the 
effect of the material conductivity on the condensing 
heat-transfer coefficient. 

5. Attach drainage strips on the existing tubes to 
Mmoceigace ene Crfece Of the drainage Strips on steam 
condensation enhancement. 

6. Take data to study the effect of vapor shear on 
finned tube enhancement. For the vacuum runs, it is recom- 
mended that vapor velocities of 1.0, 2.0, 3.0 and 4.0 m/s 
be used and for the atmospheric runs the vapor velocities 
of 0.7 and 1.0 m/s. Higher vapor velocities would be 
desirable if the apparatus could be modified. 

7. Modify the existing apparatus to study the conden- 
Sate inundation effect and take data to investigate if the 


enhancement and the optimum fin spacing are affected by 


Piunaat lon. 
8. Take data on finned tubes that are coated with a 
dropwise promoter. Since the presence of the fins will 


result in amore effective bonding of the promoter, this 

combination may result in superior enhancement and promoter 

endurance. Also, the condensate may not be retained on 

dropwise-coated finned tubes with sufficient fin spacing. 
9. Modify equation (1.14) of Owen et al. [Ref. 8] 

to include the surface-tension effect on the unflooded part 

of the tube. By including this effect, the Owen prediction 


for the steam data will become less conservative. 
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APPENDIX A 


PREPARATION OF CONDENSER TUBES 


Before every run, the tube was prepared by following 
the steps below. The operator must wear gloves to avoid 
contamination of the tube as well as to avoid skin irritation 
that may be caused by the cleaning solution. 
1. Prepare a solution of equal parts of sodium hydroxide 
and methyl alcohol. Heat the solution until it becomes 
watery. Using a toothbrush, apply a coating of this mixture 


on the Gube- 


2. Place the coated tube in a steam bath for about one 
nouiee 
3. Rinse the tube with tap water to avoid contamination 


of the apparatus with the solution. 

4. Clean the external surface of the tube with a brush 
using a weak solution of tap water and a detergent (Spark- 
leen). This will ensure that all the sodium hydroxide 
deposits are removed from the corners of the fins. 

5. Using a brush, clean the inside surface of the tube. 
Tt is essential to remove any fouling. 

6. Install the tube as soon as possible to avoid any 


system contamination from the environment. 


Pe 


APPENDIX B 


SYSTEM START-OP AND SHUT-DOWN PROCEDURES 


To start the system, the following procedure is 
recommended: 

l. Ensure that the boiler water level in the boiler is 
at least six inches above the upper level of the heaters. 

2. Close the venting valve and open the valve on the 
purge line (see Figure 2.1). 

3. Ensure that the vacuum pump is filled with oil at the 
right level. 

4. Turn on the vacuum pump. 

5. If the run is for low pressure then turn on the boiler 
power when the system pressure approximates the desired 
pressure. Avoid boiling if the system pressure is greater 
by 50 mmHg than the desired pressure. Boiling at an ele- 
vated pressure causes the apparatus walls to reach a 
temperature higher than the desired saturation temperature. 
If this happens, a considerable time interval (as much as 30 
minutes for a temperature difference of 10 K) will be needed 
to avoid superheated steam. Note that when the system walls 
are at an elevated temperature, the steam gains heat, thus 
causing superheating. The power input is controlled by a 
potentiometer until the voltmeter reads 90 units (note that 


each unit represents 2.2 Volts). This setting yields a 2.0 


in 3 


m/s vapor velocity at a system pressure of 85 mmHg. If the 
run iS an atmospheric one, then, turn on the power right 
after the vacuum pump is energized. For a vapor velocity 
of 1.0 m/s, the voltmeter should be set at its maximum 
which is 170 units” on 17030222) Onaeee 

6. Open the valve to fill the tap water sump. 

7. Circulate tap water through the auxiliary condenser. 

8. Check that the water sump is filled and then turn on 
the pumps that circulate water through the test bue. 

9. When the steam reaches the tube (as can be seen from 
the viewport), adjust the water flowrate through the auxili- 
ary condenser in order to achieve the desired pressure. 

10. Energize the data acquisition system, the computer, 
the printer and insert in the driver of the computer the 
disk which contains the code "DRP4." Finally, load "DRP4" 
and follow the interactive program to gather and analyze 
data. 

To secure the system: 


1. Turn off the input power to the boiler. 


2. Secure the water supply to the auxiliary condenser. 

3. Open the venting valve to relieve the pressure in the 
apparatus. 

4. Turn off the printer, the data acquisition system 


and the computer. 
5. Secure the pumps that circulate tap water through 


the test tube. 
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APPENDIX C 


RAW DATA OF UNINSTRUMENTED AND INSTRUMENTED TUBE 


The following two tables contain the raw data taken for 
the uninstrumented and instrumented tubes under vacuum and 


near-atmospheric conditions. 


TABLE V 
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APPENDIX D 


RoW on a Or CUBES tEoTED 


The following pages contain tables with the raw data 
taken for all tubes tested under vacuum and near-atmospheric 


Poma1tions. 
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APPENDIX E 


ERROR ANALYSIS 


The numerical data taken during an experiment are used 
together with theoretical formulations in order to obtain 
results about the objective of the experiment. Any measure- 
ment has an uncertainty which is a function of the device 
accuracy as well as the experience of the operator. 
Therefore, the final result may be severely distorted due 
to the error propagation during the computational steps. 

In cases where the final result has a large uncertainty, 
it may be unwise to draw any conclusions. 

The probable uncertainty in a computation can be 


determined from the following equation, as taken from 


[Ref. 19]: 


72 
aR 2 aR 2 JR 2 
+ (So wo) + eee + (Sw) *] 


(Eos) 


where: 


Re 1S the tess 


X,,X,,-..,X. are the measured (independent) 
variables 
W.,Wa,.--,W. are the uncertainties in the 


méasured (independent) variables. 


The above equation, together with the uncertainty of 


each measurement, was used to determine the possible error 
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of the outside heat-transfer coefficient. As discussed in 
Chapter V, the outside resistance was found by subtracting 
the water-side and tube-wall resistances from the overall 
resistance. This means that the error in the condensing 
coefficient is mainly a function of the overall and the 
water-side resistances. In such cases, if the water-side 
resistance has a significant error it will be unavoidably 
reflected in the condensing heat-transfer coefficient. 
Further, if the value of the water-side resistance is much 
larger than the condensing resistance, then the error that 
will appear in the condensing coefficient will be highly 
amplified. The steam-side resistance is generally smaller, 
especially on finned tubes, then the water-side resistance 
(recall that the condensation process is a phase change 
phenomenon where the heat-transfer process is very effec- 
tive). Therefore, to 'protect' the final result of the 
condensing coefficient from an excessive error, the water- 
Side and condensing resistances should be as nearly balanced 
as possible. The water-side resistance can be reduced 
either by increasing the coolant velocity or by internal 
enhancement. Both ways promote turbulence and reduce the 
boundary layer thickness on the inside wall surface (recall 
that the coiled insert was used exactly for this purpose). 
However, since the water-side temperature rise decreases 
with increase in water velocity, the water velocity must not 


be too excessive. Since the water-side temperature rise 1s 


al: 


the most crucial measurement, a minimum value must be 
assigned to minimize experimental errors. For this reason, 
the internal enhancement is preferable than the increase 

in water velocity to obtain a lower thermal resistance on 
the waterside. 


Based on the above discussion, it should be expected 


that as the cooling water flow rate decreases (i.e., less 
turbulence), the error in the outside heat-transfer coeffi- 
cient should increase. Figures E.1 and E.2 show the error 


bands for three cooling water flow rates for the tube that 
showed the best heat-transfer performance (1.0 mm fin thick- 
ness, 1.5 mm fin spacing and 2.0 mm fin height) for both 
pressure conditions. It is shown that indeed the error 
gets larger as the cooling water flow rate (i.e., the heat 
flux) decreases. Also notice that the error bands shown 
with the atmospheric run are smaller than the corresponding 
value for the vacuum runs. This improvement can be ex- 
plained by the increased water-side temperature rise. 

Note that the relative error decreases as the water-side 
temperature rise increases. 

Based on the foregoing discussion, it is clear why the 
reference heat fluxes for computing the enhancement ratio 
of the cross plots in Chapter V were selected as high as 
possible “033 S5randae0 MW /m* for vacuum and atmospheric 
runs, respectively). 

The uncertainty of the major variables that govern the 


condensation process (cooling water velocity, Reynolds 
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number, heat flux, LMTD, water-side and outside heat- 
transfer coefficients) are given as percentages in the 
BobovinG COMpUter OUtTpUtS Of the program "ERROR." At 
the end of this appendix, the listing of the program 


"ERROR" is provided for the Error Analysis. 


ilies ys 


DATS FUR THE ERRUR ANALYSIS: 

Powe emnes ee ee ao) 
Pressure Condition: YACUUM i? kPa 
Cooling Water lemperacuges pee 
Cooling Water lemperacinte: Hic eur 
Steam lemperature 7 12 3a0 
Nace re | en aan = $0.00 
Nel Ge tee y euler aiey, - a 
EVROR ANALYSTS: 

7, aan Dies o eee lene = {,i3 
Se On en Raae lie Number = ieee 
% EYEOr immeat Flame =e 
error aneeriee = 2,22 
, acronr Laneue ee 
/ Ertor lined = 1,3? 
Hoes ie pepe an are 2B, Sle 
DATA FOR THE ERROR ANALYSIS 

Pie theme. FeSe2Vi Seb 

7 Pessitre “Cona | Elton. ee Ole! 
COOLING Water Jenperavilze sim en 
Cooling Water Tem@eracine lt =a oeees 
Sloan Temperature =O ea 
water rlow Rate (4) =) Seay 
ote, cerows ny = oan 
SE Rie teee oe 

ee eee het ee ete eee eters = sia 
A Er¥or im Reyne! ce Numoeer ee es) 
7, Leer aimee al eee eeix = eye 
, =fror menelen = 1.63 
J, Ea fia eee et Si eee 
Jom apy vee wih 5 = eee 
jo =r mie ane =! (re, 
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DATA FOR THE ERROR ANALYSIS: 


ee Name: FeSe2A128 
mmessiire Condition: ATMUSPHER 
Cooling Water Temperature In 
Peeing Hater l[emperature Unt 
foeam ‘emperature 

memoer Flow Rate (7) 

Pewee Velocity 


ERROR ANALYSIS: 


eercor in Mass Flow Rate 
feerror in Reynolds Number 
Meeeror in Heat F lux 
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meee tror in Uo 

Peeerror in Hi 

Peer ror in Ho 
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FILE NAME : ERROR 

REVISED : JUNE 27. 1984 

PRINT USING "10X,."°DATA FOR THE ERRGR ANALYSIS:""” 
PRINI 

Beer 

INPUT “ENTER FILE NAME” ,.FilesS 

eS USING "10X.""File Name: oe olen PL les 
EEP 

INPUT “ENTER PRESSURE CONDITION ¢1=VACUUM.2=ATMOSPHERIC) °.Prce 
eae THEN 

BEEP 


PRINT USING "10X,""Pressure Condition: VACUUM 11 kPa?” 

PRINT USING “°10X,""Pressure Condition: ATMOSPHERIC 101 kPa’’” 
END IF 

BEEP 

INPUT "ENTER COOLING WATER TEMPERATURE IN". Tci 


anievolNG  1T0X," “Cooling Water Temperature Im o= “" ,30.20,°° (Deg ©)": 
Beer 

INPUT “ENTER COOLING WATER TEMPERATURE DIFFERENCE” ,Dtw 

Tco=lIcitDtu 

See USING 10%. Cooling Water Temperature Qut = ~*,.30:2), °° (Deg CG)”; 
BEEP 

INPUT “ENTER STEAM TEMPERATURE’’.Ts 

PRINT USING "“10X,""Steam Temperature ee Oe D 9 Vogme) 
BEEP 

INPUT "ENTER COOLING WATER FLOW RATE Z%*',Fm 

meeuNi USING “10X,." "Water Flow Rate (4) = oun AF im 


T=(TcitTco)/2 ' FILM TEMPERATURE ° 
UNCERTAINTY IN THE COOLING WATER 

Dfm=.10 ' ERROR IN FLOW METER READING 

Drho=1 ! ERROR IN WATER DENSITY 
Dmf=5.00049E-3+6.3861937E-3*Dfm *' ERROR IN MASS FLOW RATE 
Rho=FNRho(T) ! WATER DENSITY 

Mf=1 0480SE-2+6.30932E-3*Fm * MASS FLOW RATE OF COOLING WATER 
CORRECT MF FOR THE TEMPERATURE EFFECT 

MF =MFx¢(1,.0365-1.96644E-3*TcitS. 252E-6*Tc1 °2)/.995434 

Di=.0127 ¢' TUBE INSIDE DIAMETER 

Ai=(PI*Di~2)/4 ! TUBE INSIDE CROSS SECTION AREA 

Ddi=.000025 
Dar=(PI*Ddi°2)/4 ' ERROR OF INSIDE TUBE CROSS AREA 
COMPUTE THE WATER VELOCITY 

Vu=MFf/CRho*Ai) |! WATER VELOCITY 

PRINT USING "10X,""Water Velocity = ae. oe MMS) ot 
CORRECT OUTLET WATER TEMP. FOR THE MIXING CHAMBER EFFECT 

Tco=Tco-.004*Vw°2 '| COOLING: WATER OUTLET TEMP. 

imei t+ico)*.5 ! FILM TEMPERATURE 

COMPUTE THE ERROR IN WATER VELOCITY 

Dvw=Vu*SQR¢( (Dmf /MF) >2+¢(Drho/Rho) -2+(Dai/Ai)*2) 

UNCERTAINTY IN THE REYNOLDS NUMBER 

Mw=FNMw(T) § WATER VISCOSITY 

Dmw=6.E-6 '! ERROR OF WATER VISCOSITY 

Re=(Rho*Vw*Di >) /Miy 
Dre=Re*SQR((Drno/Rho) 2+ (Dvw/Vw) °2+¢Ddi/D1)°2+¢Dmu/Mw) ~2) 


roo 
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1560 
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1580 


UNCERTAINTY IN 
Cow=FNCpuw(T) 
Q=Mf«(Tco-[c1) Cow 

Dtco=.U2 

Dter=.02 

Dcpw=3 

Da=O*SQR( (Dmf /MF) > 24+¢¢Dtco/( Teo-Te:1))) “2+¢(Dte1/(Tco-Te1)))°2+¢(Depu/Couw) ~2 


THE BEAT TRANGSE Srey 


* UNCERTAINTY IN THE HEAT FLUX 


Dl=.0005 ! ERROR IN TUBE LENGHT 

Do=.01905 *' SMUOTH TUBE OUTER DIAMETER 
Ddo=.000025 

L=.13335 ! FINNED TUBE LENGTH 

Qp=Q/¢(PI*Do*xL) *' HEAT FLUX 
Dap=(Qp/PI)*SQR¢(Dq/Q)°2+(Ddo/Do)*2+(DIL/L) *2) 
Dts=.02 !' ERROR IN TEMPERATURE MEASUREMENT 
Lmtd=(Tco-Tci1)/LOG¢¢Ts-Te1)/¢(Ts-Teco)) 
Uo=Qp/Lmtd ' QVERALL HEAT TRANSFER COEF. 
Alt=Dts*(Tci1-Tco)/((Ts-Tei)*¢(Ts-Tco)#LOG(¢(Ts-Te1)/¢Ts-Teo) >? 
A2=Dtc1/¢(Ts-Te1)*LOG((Ts-Tci)/(Ts-Jco))) 
A3=Dtco/((Ts-Tceo) *LOG¢(¢(Ts-Te1)/¢(Ts-Tco))) 
Dimtd=Lmtd*SQR¢(A1*2+A2°2+A3° 2) 
Duo=Uo*SQR¢((Dap/Qp)°2+¢(Dilmtd/Lmtd) ~2) 

M=My 

P1201 4276815972732. /s 

Kw=FNKw(T 1) 

Ci=.0635 ! SIEDER-TATE COEFFICIENT 

Ac=26.4 ? INTERSCEPT FROM SIEDER PROGRAM 
Kce=385 ' WALL THERMAL CONDUCTIVITY 

L1=.060325 ! LENGTH OF UNFINNED LEFT PART OF TUBE 
L2=.034925 ! LENGTH OF UNFINNED RIGHT PART OF TUBE 
Dr=.015875 

Pr=Cow=Muy/Kuy 

Muw=F NMuw¢ T) 

UNCERTAINTY OF INSIDE HEAT-TRANSFER COEFF. 
Cf=1. , 
Hi=(Kw/Did*(CixRe .8*Pr>.333*Cf+Ac) 
Dti=Q/¢(PI*Dix(L+L1*Fel+L2*Fe2) *H1) 
Cfco=(Muw/FNMuw(T+Dt1))>.14 

IF ABS((CFCE-CF)/CFCc)>.01 THEN 

Cf=(Cf+Cfe)*.5 

GOTO 1900 

END IF 

P1=PI*(Dit+Do) 

B1=(Do-Di)*PI=~(Di+Do)=.5 

M1=(Hi¥P1/(Kco~81))°.5 

P2=PJT=(DitDr>) 

B2=(Dr-D1)«PIl*(Di+Dr)~=.5 

M2=(HixP2/ (Kce*B2))°.5 

Fel=FNTanh(M1*L1)/¢(M1«L1) 
Fe2=FNTanh(M2*L2)/(N2*L2) 
Dtc=Q/(PI*Dix¢(Lt+L1*Fel+L2*Fe2) «Hi ) 

[IF ABS((Dte-Dt1)/Dtce)>.91 THEN 1900 

Dkw=.0010 ! ERROR IN WATER THERMAL CONDUCTIVITY 
Dei=.0005 ! ERROR IN SIEDER-TATE COEFFICIENT 
Dpr=.05 * ERROR IN PRANDTL NUMBER 

Def =8.E-6 

A4=.14*Dcf/CFf 

Dhi=HixSOR¢((Dkw/Kw) ~“2+(Ddi/Di)d 2+¢(.2@=*Dre/Re) “2+ .333*Dpr/Pr) “2+(De1/C1) “2+ 


UNCERTAINTY OF OUTSIDE HEAT-TRANSFER COEFF. 
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DOwr~ NN SE 
OoOOwr DOO @® 


meeeccULATE THE “£ ERROR 
f CALCULATE THE % ERROR 
PecALCULATE THE 4 ERROR 
meenlCULATE THE Z ERROR 
meeneGULATE THE “7 ERRGR 


Dkc=1 ' ERROR IN TUBE THERMAL CONDUCTIVITY 

Ru=Do*_UG(Do/D1)/(2*Ke) §' WALL RESISTANCE 

Ho=17°(¢1/U0)—(Do#Ll/(D1*(L+t+Li*fel +L2*Fe2)*Hi)d)-Rwy) 

Drw=Rw*SQR( (Ddo/Do) °2+¢(Dkce/Ke) “2+(Ddo/ (Do#*LOG(Do/D1)))°2+¢(Ddi1/(D1*LUG(Do/D 


“2) 


Ao=!/U0o-Rw-(Do/ (D1 *Hi)>) 
A6=Duo/ (U0 2*AS) 


A/7=Drw/45 
A8=((Do/(DixHi) )*¢(Dhi/Hi))/AS 
PRINT 
Dho=Ho*SQR(A6B72+A772+AB*2) 
f CALCULATE THE Z% ERROR IN Ho 
Prho=Dho*#100/Ho 
' CALCULATE THE Z% ERROR IN REYNOLDS NUMBER 


Prre=Dre*100/Re 
IN 


IN 
IN 
IN 
IN 


MASS FLOW RATE 

HEAT TRANSFER 

ere 

QVERALL HEAT TRANSFER COEF. 
PSE eG SIkGNSrER CHER 


Prmt=Dmt «100/MFf 
Prap=Dqp*100/Qp 
Primtd=Dlmtd*100/Lmtd 
Pruo=Duo*100/Uo0 
eum -On1*100/H1 


PRINT 

PRINT USING "10X,"”"”ERRGR ANALYSIS:'**** 

PRINT 

PRINT USING ‘10X,"°°% Error in Mass Flow Rate =e 2D. eD et erm 
PRINT USING "10X,"°"°% Error in Reynolds Number =oencln cower Te 

PRINT ISSING °10X."°°% Error in Heat Flux 2 7 De P Bap 
PRINT USING "10X."°"°% Error in LMTD = "2.20" :Primtd 

PRINT USING “10X,"°"°% Error in Uo =e en ee LUO 
PRINT USING "10X,.""°% Error in Hi =— ee cD, ETAL 

met USING "10X,"°"°% Error in Ho = Feviecl. sero 
N 


DEF FNMuw¢T) 

A=247 .8/¢(7T+133.15) 
Muw=2.4E-5*10°A 

RETURN Muw 

FNEND 

DEF FNTanh(X) 

P=EXP(X) 

Q=EXP¢(-X) 

Tanh=(P-Q)/(P+Q) 

RETURN Tanh 

FNEND 

DEF FNKuw(T1) 
Kw=-.92247+7T1#«(2.8395-T1«¢( 1 .8007-T1#*¢.52577-.07344*T1))) 
RETURN Ky 

FNEND 

DEF FNMu(T) 

B=747 .8/¢(1T+133.15) 
Mw=2.4E-5*10°A 

RETURN Mw 

FNEND 

DEF FNRho(T) 

Rho=999 .52946+T*¢ .01269-T*(5.482512E-3-T#!.234147E-5)) 
RETURN Rho 

FNEND 


eta 


2/20) Dem E NEP CT 
2/40 Cow=(4,21120858-1<(2. 26826E-3-12(4, 42361-5342) 42 ce 7 een 


2750 RETURN Cpw 
2760 FNEND 
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